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A COMPREHENSIVE COMPUTER PROGRAM FOR PREDICTING 


SOLAR CELL PERFORMANCE 
J. P. Kirkpatrick, C. A. Lewis, and J. C. Caspar 
Ames Research Center 


SUMMARY 


A comprehensive computer program has been developed to predict the performance of 
individual solar cells and complete arrays as a function of solar intensity. The program is unique in 
its capability to calculate the light-generated current of a single cell by integrating the solar or 
simulator spectrum with the following parameters as a function of wavelength: solar cell spectral 
response, filter transmission, and space radiation damage. (Variations with temperature and 
incidence angle can also be handled.) The diode equation is then used to obtain the current-voltage 
characteristics and maximum power point. Using a variable curve factor derived from experimental 
data in the diode equation enhances the prediction accuracy over a wide range of light intensities. 
Array current-voltage curves may be calculated for combinations of series and parallel cell wiring 
arrangements. These curves reflect wiring and diode losses as well as angular orientation to the sun. 

The program is also useful in extrapolating available data, comparing experimental data taken 
with different energy sources, and correlating experimental and analytical predictions. 


INTRODUCTION 


Optimization of solar array design and performance prediction for future missions is becoming 
more complex as spacecraft encounter wider variations in the solar environment and extended 
exposures to damaging space radiation. The use of thermal control techniques to provide 
satisfactory cell temperatures may also introduce variables that must be considered. For example, a 
parametric study of solar array performance for a near solar mission may require that 
coverglass-filter transmission, solar cell response, and space radiation damage be considered on a 
wavelength basis as a function of solar distance, angle of solar incidence, and cell temperature. 
Correlation of data obtained with different energy sources also requires spectral calculations. 

The program described, unlike current programs (refs. 1-4), is capable of predicting cell 
performance over a wide range of light intensities using spectral information directly. It is written in 
Fortran and is operational on both the IBM 7040/7094 Direct Couple System (DCS) and the IBM 
360/67 computers with optional graphic output on an S-C 4020. 

To provide maximum utility, the program calculates individual cell or total array performance 
on the basis of either a few basic cell parameters or detailed spectral data which may be a function 
of one, two, or three independent variables. A modification of the classical diode equation (ref.. 5) is 
used to construct the current-voltage curve from which maximum power is determined. A variable 
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curve factor derived from experimental data is used to enhance the accuracy of the diode equation 
in predicting the shape of the current-voltage curve at high hght intensities. 

Total array performance is synthesized by combining single cells into parallel and series wiring 
configurations to form strings and then accounting for the total number and solar orientation of 
each string in the array. 


SOLAR CELL MODEL 


Current Density 

Incident light energy is converted by a solar cell into electrical current. The short circuit 
current density is the current flow per unit of cell area developed at the terminals when the external 
load has zero resistance. It is determined from the following expression: 

JCD = 140 cos THETA {xjlTR[1 + ATX(T - TREF)]} (1) 

AU2 

where 

JCD short circuit current density, mA/cm* 

THETA angle between the solar vector and the normal to the cell surface, deg 
AU solar distance, AU 

XJLTR conversion constant at reference conditions, mA/mW of incident energy 

ATX temperature coefficient defined as the dimensionless change in current per °K change 

in temperature. When XJLTR is calculated by method 2 (discussed below), ATX is 
for the bare ceU when filtered coverglasses are being considered. 

T cell operating temperature, °K 

TREF reference temperature for standard test conditions, typically taken as 301 ® K 

Equation ( 1 ) states that the maximum current density is determined by the product of the normal 
component of the incident light energy and the power-to-current conversion constant as modified 
by the deviation in temperature from reference conditions. 

The temperature coefficient is usually provided by the cell manufacturer either directly or in 
curve form. The conversion constant, however, may not be available, and under such circumstances 
the computer program will calculate this factor by one of two different methods. The first requires 
a measurement of short circuit current (IREF) under standard text conditions (i.e., light 
intensity = 140 mW/crn^ at temperature TREF). This is usually available from the manufacturer. 
The computer solves the equation 
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XJLTR = 


IREF 


( 2 ) 


AREA X 140 

where AREA is that of the solar cell. XJLTR is then considered a constant that does not change as 
parameters deviate from standard test conditions. At large angles of incidence and temperature 
extremes, this assumption may be invalid as a result of changes in coverglass-filter transmission and 
cell spectral response. The second method is suggested for calculating XJLTR under these 
conditions. 

The second method constitutes a unique feature of this computer program and evaluates 
XJLTR by integrating on a wavelength basis the incident energy, the solar cell response, and the 
filter-coverglass transmission (if used). Both the transmission and response characteristics may be 
degraded for space radiation damage. This integration is performed by the following operation: 


XJLTR = 


1.3 

CONST ^ {Lj^[Rj^ (THETA, T) - ARj^ (AU, THETA) ] [TRj^ (THETA, T) - ATR^ (AU, THETA)] } 
0.3 


(3) 


where 


Ljy^ relative spectral energy distribution of the light source expressed as a ratio of energy in 

(LSOLAR) each wavelength interval (usually 0.01 n) to the total energy (This information is 
interpolated from tables LSOLAR(L\) and W A VE( wavelength), defined by data 
statements in subroutine RESPE. These tables are of length NW. The wavelength 
values in this statement are used as the discrete summation points in the numerical 
integration.) 


R\ relative spectral response of the solar cell in dimensionless units in the form of an 

(RESP) interpolated table with values of RESP as a function of wavelength; or wavelength 

and incidence angle; or wavelength, incidence angle, and temperature 


TR\ spectral transmission expressed as a ratio of transmitted to incident light for the 

(TRANS) coverglass-filter assembly, in the form of an interpolated table with values of 
TRANS as a function of wavelength; wavelength and incidence angle; or wavelength, 
incidence angle, and temperature (For bare cells, this table may be eliminated and 
TRANS is automatically set equal to unity.) 


AR\ degradation in cell response caused by space radiation from the undegraded values 

(DRESP) listed in the response tables, in the form of an interpolated table with values of 
DRESP as a function of wavelength; wavelength and incidence angle; or wavelength, 
incidence angle, and AU (position in orbit) 

ATR\ degradation in coverglass-filter transmission caused by space radiation from the 
(DTRANS) undegraded values hsted in the transmission tables, in the form of an interpolated 
table with values of DTRANS as a function of wavelength; wavelength and incidence 
angle; or wavelength, incidence angle, and AU (position in orbit) 
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CONST proportionality factor that converts the relative cell response used by the numerical 
integration to absolute units of mA/mW. Its value may be input directly or 
calculated by the program from the following equation: 

CONST = (4) 

AREA X TIRR C^REF;^) (R;^) ] 

X 

where 

IREF short circuit current in mA at reference conditions (TREF, LREF, TIRR). This is 
normally at air mass zero (amz). 

LREF\ relative spectral energy distribution of the source used to obtain IREF (Johnson’s solar 
(LREF) spectrum (ref. 6) is used by the program, but is easily replaced if the 

associated amz short circuit current is not available. LREF^ is defined in 
subroutine RESPE by data statement LREF, the values of which correspond to 
center wavelengths defined by data statement RWAVE. The tables are of length NP.) 

AREA cell area, cm^ 

TIRR total irradiance of source used to obtain IREF in mW/cm^ 

R\ relative spectral response of solar cell at TREF and zero angle of incidence 

(RESP) 

Once XJLTR has been determined and equation (1) solved, the light-generated current (IL) is 
obtained by the product of the current density (JCD) and the cell area, that is, 

IL = JCD X AREA 

XJLTR is the dominant parameter affecting the value of JCD and an uncertainty in its 
magnitude should be reflected in IL. ATX, however, is generally less than 0.001 and becomes 
important only for large deviations from the reference temperature. 

Open Circuit Voltage 

The open circuit voltage (VOC) is determined from the following expression: 

VOC = VREF - BRS(T - TREF) + VCL(IL - IREF) (5) 

where 

VREF open circuit voltage under standard test conditions, mV 

BRS voltage-temperature coefficient defined as the change in voltage per degree change in 
temperature, mV/°K 
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VGL voltage-light intensity coefficient defined as the change in voltage per change in 
light-generated current from standard reference conditions, mV/mA 

T, TREF, IL, and IREF are as previously defined. 

The voltage-temperature coefficient (BRS) is an important parameter and is generally available 
from the manufacturer. The voltage-light intensity coefficient (VCL) is derived from experimental 
data and becomes significant only for large variations in light intensity. Typical values of VCL are 
included in the sample problems. 


Basic Diode Equation 

The program uses the classical diode equation for a solar cell to generate current -voltage (I-V) 
curves for specific operating conditions. This equation may be written as; 

1 . 10 Lxprqf'' - I ■ I^S)] . 

I *^L 1000 AxKxT J 

where 

I load current, mA 

10 reverse saturation current, mA 
V load voltage, mV 

RS series resistance, 

IL light-generated current, mA 

A factor used to fit experimental curve to theoretical curve 
K Boltzmann’s constant, 1.38X1 J/°K 

T temperature, °K 

q electronic charge, 1.60X1(T*® coulombs 

The above expression is based on the generally applied, single-stage equivalent circuit for a 
solar cell. To improve prediction accuracy over a wide range of light intensities, a variable curve 
factor A derived from experimental data is employed. The optimum A factor as a function of 
light intensity is depicted in figure 1 and has been found to yield accuracies within 3 percent at the 
maximum power point for cells of various manufacturers tested at solar intensities up to 25 solar 
constants. If accurate predictions are required at high light intensities without a variable curve 
factor, it would be necessary to use a more complicated model for the solar cell utilizing distributed 
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parameters for a more realistic representation of the circuit operation. This, in turn, would result in 
a complex mathematical expression for the current-voltage relationship requiring basic cell 
information not readily available. 

The series resistance (RS) should be available from the cell manufacturer, but may be 
computed by the measurement technique of Wolf and Rauschenbach (ref. 5) if necessary. 

The standard reverse current for a diode; 10, is determined from a solution of the diode 
equation under open circuit conditions. This gives the following: 


10 


exp 


( 


IL 

qVOC 

1000 AxKxT 


) 


(7) 


The computer program first solves equation (5) to obtain VOC. Then 10 is determined by 
equation (7). Incremental values of V are then used in equation (6) to generate an I-V curve. The 
maximum power point is then computed from these data. 


ARRAY CALCULATIONS 


Performance may be calculated for solar arrays that are either sun-oriented or spinning and 
composed of a single flat paddle or a complex combination of panels with different angular 
orientations to the sun. 

The string current-voltage (I-V) characteristic is first determined at each desired angular 
orientation (THETA) by multiplying the current coordinate (for a series of incremental voltages) of 
the individual cell I-V curve by the number of cells connected in parallel (NIS) and the voltage 
coordinate by the number of cells in series (VIS). The voltage coordinate thus obtained is then 
reduced by the string circuit loss (SCL). Algebraically, this may be represented by 


^(^^string NIS X I(i)cell 

V(i)string = VISX V(i)ceii-SCL (9) 


where (i) varies in steps from zero to the single-cell open circuit voltage (VOC). 

The current coordinates of the string I-V curves are next multiplied by the number of strings 
(FTHETA) at each angular orientation of the sun, and the current contributions at each angle are 
then summed to obtain the final shape of the array 1-V curve. String and array voltage coordinates 
are assumed to be identical. These operations may be represented by 
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I(i) 


array 


= string ■< ™™] 


THETA 


V(i) = V(i) ^ . 

array string 


( 10 ) 

( 11 ) 


where (i) varies as before. 

The maximum power output of the solar array is then determined from the array I-V curve. 


PROGRAM FLEXIBILITY 


The program accepts input information ranging from single-valued coefficients to complex 
trivariant tables. It can predict the performance of a single cell or a multifaced array throughout an 
entire mission. Characteristic of the varied tasks performed by the program are the following 
calculations that appear in detail as sample problems in appendix A. 

1. Single-cell current-voltage curve and maximum power point determined from a minimum of 
cell parameters with graphical output of power versus voltage and current versus voltage 

2. Array I-V characteristics and maximum power point for a cylindrical array with its spin 
axis normal to the ecliptic plane 

3. Filtered cell output calculated from coverglass-filter transmission and cell response input as 
a function of wavelength 

4. Performance of a filtered cell during laboratory tests using simulated solar radiation 

5. Complete array calculations, including space radiation damage, for a spinning spacecraft 
with solar panels inclined at an angle to the solar vector. The cell response is a function of 
wavelength and incidence angle, and the coverglass-filter transmission is a function of wavelength, 
incidence angle, and temperature. 

In addition to the mission-oriented applications described above, the program is also well 
suited to basic studies of various cell parameters such as series resistance, temperature coefficients, 
and spectral response. 


CONCLUDING REMARKS 


Basic cell and total array performance can be predicted for a wide range of environmental 
conditions using known coefficients and either test data taken under standard reference conditions 
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or spectral information on cell response, light source, and (if desired) coverglass-filter transmission 
and space radiation damage. The accuracy of computed I-V curves and maximum power points is 
commensurate with the accuracy of the input data. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., 94035, Apr. 1, 1970 
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APPENDIX A 


SAMPLE PROBLEMS 

SAMPLE PROBLEM 1 - SINGLE-CELL CALCULATION 


A particular low-series resistance solar cell is being considered for a Mercury mission. The 
output characteristics of this cell at 0.50 AU, normal incidence, and a panel temperature of 360° K 
are desired. The following data are available from the manufacturer’s information and laboratory 
tests: 

IREF 60 mA (short circuit current at 1 solar constant) 

VREF 594 mV (Open circuit voltage at 1 solar constant) 

TREF 303° K 

RS 0.11 (cell series resistance) 


light intensity, 
sun 

Temperature (T), 
°K 

Current (I), 
mA 

Voltage (V), 
mV 

1 

303 

60.0 

594 

1 

393 

61.8 

397 

6 

393 

376 

464 


The coefficients ATX, BRS, and VCL may be calculated from these data as follows: 


ATX 


~ = = 0.00035/°K 

T 2 - Tj 60.0 393 - 303 


BRS 


Vj - V 2 _ 594 - 397 
T 2 - Ti " 393 - 303 


2.0 mV/°K 


VCL 



464.0 - 397.0 

376.0 - 61.8 


= 0.21 mV/mA 


Figure 1 indicates a value of 1 .43 for A. Detailed output and plotting are desired, but radiation 
damage is not to be considered. 


The input data form, data deck listing, and computer output are shown in figures 2, 3, and 4; 
current-voltage and power-voltage curves are plotted in figures 5 and 6. 
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Figure 1.- Variation in "A" factor with intensity. 
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INPUT DATA - SOLAR CELL PERFORMANCE PROGRAM 


FORMAT INPUT DATA 
A6 RUN I 

I2A6 TITLE CARD Sample Problem I - Single-Cell Calculation 

8110 NAU 1 NTHETA 1 KPLOT 1 KOUT 1 INT 0 INR 0 INTD 0 INRD_0 

8F10.4 AU 0.50 

TEMP 560 ■ 

A 1.43 

VCL 0.21 

AREA 2 ^ TREF 303. VREF 594. I REF 60.0 XJLTR Q. ATX .00035 
BRS 2.0 RS 0.11 

THETA 0 . 


FTHETA 1 . 

NIS 0 . VIS SCL 

PLOT CARD 08/30/68 KIRKPA AMES 244 

TABLES 

TRANSMISSION (TRANS) Not used - delete 


RESPONSE (RESP) Not used - delete 


TRANSMISSION DAMAGE (DTRANS) Not used - delete 


RESPONSE DAMAGE (DRESP) Not used - delete 


COMMENTS Single-cell calculation with a minimum of required input 
information 


Figure 2.- Input data form - Sample Problem 1 
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RUN I 


SAMPLE PRObLEM 1 - SINGLE CELL CALCULATION 
111 100 

.50 

360. 

1.^3 

0.21 

2.0 303. 594. 60.0 0. .00035 

0 . 

1 . 

0 . 

09/16/69 J CASPAR AMES 233 


Figure 3.- Data deck listing - Sample Problem 1 


0 0 


2. .11 


INPUT DATA RUN 1 


SAMPLE PROBLEM 1 - SINGLE CELL CALCULATION 


NAU = 

1 


NTHETA= 

1 


KPLOT= 

1 


KUUT = 

1 


AREA = 

0.20000000E 

01 

TKEF = 

O.3O3OU0O0E 

03 

VREF = 

0.59400000E 

03 

IREF = 

0.60000000E 

02 

XJLTR= 

0.0 


ATX = 

0.34999987E- 

-03 

BRS = 

0.20000000E 

01 

RS = 

0.10999995E 

00 


Figure 4.- Output 


Sample Problem 1. 


1^/U1/6V 


SULAR CELL PERFURmANCE 
RUN I 


SAhPLE PRUbLEh 1 - SINGLE CELL CALCULATION 
FIXED CELL CHARACTERISTICS 


AREA* 2.00 VKEF=s 59A,0 IREF = 60.00 ATX= 0.000350 bRS= 2.00 RS* 0.11 


ENVIRONMENTAL CELL PERFORMANCE - (SINGLE CELL) 


AU= 0.50 theta = 0.0 TC =360.0 JLTR = 0.21^28567E 00 VCL =0.21 
IL ^ 0.2A478790E 03 10 = 0. 2062 1 1 19E-02 VUC = 0.51880542B 03 A = 1.43 


I 

CURRENT (MA ) 

VULTAGE(MV) 


1 

0.24478790E 

03 

0.0 


2 

0.24478790E 

03 

O.IOOOOOOOE 

02 

3 

0.24478790E 

03 

0.20000000E 

02 

4 

0.24478790E 

03 

0.30000000E 

02 

5 

0.24478790E 

03 

0.40000000E 

02 

6 

0.24478790E 

03 

0.50000000E 

02 

7 

0.24477534E 

03 

0.60000000E 

02 

8 

0.24477534E 

03 

0.70000000E 

02 

9 

0.24477534E 

03 

0.80000UOOE 

02 

10 

0.24476125E 

03 

0.90000000E 

02 

11 

0.24475401E 

03 

0. lOOOOOOOE 

03 

12 

0.24475314E 

03 

O.llOOOOOOE 

03 

13 

0.24473355E 

03 

0. 12000000E 

03 

14 

0.24471930E 

03 

0.13000000E 

03 

15 

0.24471017E 

03 

0. 14000000E 

03 

16 

0.24467911E 

03 

0.15000000E 

03 

17 

0.24465112E 

03, 

0. 16000000E 

03 

IB 

0.24462579E 

03 

Oa7000000E 

03 

19 

0.24457216E 

03 

0.18000000E 

03 

20 

0.24451720E 

03 

0.19000000E 

03 

21 

0.24444833E 

03 

0.20000000E 

03 

22 

0.24436487E 

03 

0.21000000h 

03 

23 

0.24425423E 

03 

0.22000000E 

03 

24 

0.24411565E 

03 

0.23000000E 

03 

'2 5 

0.24394846E 

03 

0.24000000E 

03 

26 

0.24373842E 

03 

0.25000000E 

03 

27 

0.24347366E 

03 

0.26000000E 

03 

28 

0.24314253E 

03 

0.27000000E 

03 

29 

0.24272852E 

03 

0.28000000E 

03 

30 

0.24221111E 

03 

0.29000000E 

03 

31 

0.24156488E 

03 

0.30000000E 

03 

32 

0. 240758 30E 

03 

0.31000000E 

03 

33 

0.23975249E 

03 

0.32000000E 

03 

34 

0.23849962E 

03 

0.33000000E 

03 

35 

0.23694106E 

03 

0.34000000E 

03 

36 

0.23500552E 

03 

0.35000000E 

03 

37 

0.23260670E 

03 

0.36000000E 

03 

38 

0.22964098E 

03 

0.37000000E 

03 

39 

0.22598561E 

03 

0.38000000E 

03 

40 

0.22149620E 

03 

0.39000000E 

03 

41 

0.21600662E 

03 

0.40000000E 

03 

42 

0.20932799E 

03 

0.41000000E 

03 

43 

0.20125087E 

03 

0.42000000E 

03 

44 

0.19154869E 

03 

0.43000000E 

03 

45 

0. 17998454E 

03 

0.44000000E 

03 

46 

0. 16631820E 

03 

0.45000000E 

03 

47 

0.15031766E 

03 

0.46000000E 

03 

48 

0.13176140E 

03 

0.47000000E 

03 

49 

0.11048154E 

03 

0.48000000E 

03 

50 

0.86327515E 

02 

0.49000000E 

03 

51 

0.59201401E 

02 

0.50000000E 

U3 

52 

0.29058151E 

02 

0.51000000E 

03 

53 

0.0 


0.51880542E 

03 


I(PMAX) s 0.21871051E 03 V(PhAX) = 0.39534399E 03 PlhAX) = 0.86467926E 02 

Figure 4.- Concluded. 


13 




Figure 5.- Current -volt age plot (S-C 4020). 



Figure 6.- Power-voltage plot (S-C 4020). 
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SAMPLE PROBLEM 2 - ARRAY CALCULATION 


The electrical output of a solar array operating at 0.5 AU and 363° K is of interest. The cells 
are similar to the one described in Sample Problem 1 and their characteristics are assumed to be 
independent of temperature and solar incidence angle. The array is a spinning cylinder with its spin 
axis normal to the solar vector as shown in figure 7. For convenience, the panels may be considered 
to be a 36-sided polygon with surface normals at angles with the solar vector of 5°, 15°, 25°, etc. 
Only half the array is exposed to the sun at any given moment and this half will be symmetrical 
about the solar vector. 

There will be a total of 48 strings, each composed of 4 cells connected in parallel and 54 of 
these parallel groups connected in series. The circuit loss for each string is estimated to be 0.5 V. 
For each segment, there will be 48/36 = 1 .333 strings. Twice this number of active strings may be 
considered at each angle from 5° to 85° because of quadrant symmetry. 

The cell parameters were calculated in the previous problem and will not be repeated here. 
Detailed output is desired. 

The input data form, data deck listing, and computer output are shown in figures 8, 9, and 10. 
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\0 

t- 


(X) 



\6 



INPUT DATA - SOLAR CELL PERFORMANCE PROGRAM 


FORMAT INPUT DATA 
A6 RUN 2 


6110 NAU_1 

NTHETA 9 KPLOT 0 KOUT 1 

INT 0 

INR 0 

INTO 0 

INRD 0 

8F10.4 AU 

0.50 





TEMP 

360. 





A 

1.43 





VCL 

0.21 





AREA 

2.0 TREF 303. VREF 594. I REF 

60.0 

XJLTR 0 , 

ATX 

.00035 


BRS 2.0 RS 0.11 





THETA 

5. 15. 25. 35. 

45. 

55. 

65. 

75. 


85. 





FTHETA 

2.667 2.667 2.667 2.667 

2.667 

2.667 

2.667 

2.667 

2.667 

NIS 

4. VIS 54. SCL 0.5 





PLOT CARD_ 

Not used - delete 






TABLES 

TRANSMISSION (TRANS) Not \ised - delete 


RESPONSE (RESP) Not used - delete 


TRANSMISSION DAMAGE (DTRANS) Not used - delete 


RESPONSE DAMAGE (DRESP) Not used - delete 


COMMENTS Array calculation with a minimum of required input information 


Figure 8.- Input data form - Sample Problem 2 
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KUN 2 


SAMPLE PROBLEM 2 “ ARRAY CALCULATION 


1 

0.50 

360. 

1.43 

0.21 

9 

0 

1 

0 

0 

0 

0 

2.0 

303. 

594. 

60.0 

0 . 

.00035 

2. 

.11 

5. 

85. 

15. 

25. 

35. 

45. 

55. 

65. 

75. 

2.667 

2.667 

4. 

2.667 

54. 

2.667 

0.5 

2.667 

2.667 

2.667 

2.667 

2.667 


Figure 9.- Data deck - Sample Problem 2. 


INPUT DATA RUN 2 


SAMPLE PROttLEM 2 - ARRAY CALCULATION 


NAU = 

1 


NTHETA= 

9 


KPLOT= 

0 


KOUT = 

1 


AREA = 

0.20000000E 

01 

TREF = 

0.30300000E 

03 

VREF = 

0.59400000E 

03 

IREF = 

0.60000UUOE 

02 

XJLTR= 

0.0 


ATX = 

0.34999987E-03 

BRS = 

0.20000000E 

01 

RS- 

0. 10999995E 

00 


Figure 10.- Output - Sample Problem 2. 



12/01/69 


SOLAR CELL PeRh'URMANO: 
RUN 2 


sample problem 2 - ARRAV CALCULATrON 
FIXED CELL characteristics 


areas 2.00 VKEFs 59A,0 IREFsbO.OO ATX= 0.0003SO BRS= 2.00 RS« 0,11 


environmental cell performance - (SINGLE CELL) 


AU= 0»50 theta = 5.00 TC =360.0 JLTR = 0.21A28567E 00 VCL =0.21 
IL = 0.2A385635E 03 10 = 0 , 20633^3 lE-02 VUC = 0.518609626 03 A = 1,^3 


I CURRENT (MA) VOLTAGE (MV) 


1 

0,243856356 

03 

0.0 


2 

0.24385635E 

03 

O.IOOOOOOOE 

02 

3 

0.24385635E 

03 

0.20000000E 

02 

4 

0.24385635E 

03 

0.30000000E 

02 

5 

0.243856356 

03 

0.400000006 

02 

6 

0.24385635E 

03 

0.50000000E 

02 

7 

0.24384383B 

03 

0.60UOOOOOE 

02 

8 

0.24384383E 

03 

0.70000000E 

02 

9 

0.24384383E 

03 

0.800000006 

02 

10 

0.24382976E 

03 

0.900000006 

02 

11 

0.24382253E 

03 

0. lOOOOOOOE 

03 

12 

0.24382150E 

03 

O.llOOOOOOE 

03 

13 

0.24380209E 

03 

0. 120000006 

03 

14 

0.24378787E 

03 

0.130000006 

03 

15 

0.24377814E 

03 

0.140000006 

03 

16 

0.24374774E 

03 

0.15000000E 

03 

17 

0.24371982E 

03 

0. 16000000E 

03 

18 

0.243694156 

03 

0.17000000E 

03 

19 

0.24364099E 

03 

o.ieooooooe 

03 

20 

0.24358611E 

03 

0.19000000E 

03 

21 

0.24351736E 

03 

0.200000006 

03 

22 

0.24343398E 

03 

0.210000006 

03 

23 

0.24332359E 

03 

0.220000006 

03 

24 

0.24318596E 

03 

0.23000000E 

03 

25 

0.24302032E 

03 

0.24000000B 

03 

26 

0.24280867E 

03 

0.250000006 

03 

27 

0.24255042E 

03 

0.26000000E 

03 

28 

0.24221378E 

03 

O.27OOOO0OE 

03 

29 

0.24179807E 

03 

O.28OO0000E 

03 

30 

0.24128394E 

03 

0.290000006 

03 

31 

0.24063879E 

03 

0.300000006 

03 

32 

0.2'3983357E 

03 

0.310000006 

03 

33 

0.23882944E 

03 

0.32000000E 

03 

34 

0.23757864E 

03 

0.33000000E 

03 

35 

0.23602267E 

03 

0.34000000E 

03 

36 

0.23409032E 

03 

0.350000006 

03 

37 

0.23169537E 

03 

0.360000006 

03 

38 

0.22873441E 

03 

0.37000000E 

03 

39 

0.22508467E 

03 

0.380000006 

03 

40 

0.22060242E 

03 

0.39000000E 

03 

41 

0.21512106E 

03 

0.40000000E 

03 

42 

0.20845224E 

03 

0.41000000E 

03 

43 

0.20038654E 

03 

0.420000006 

03 

44 

0.19069771E 

03 

0.430000006 

03 

45 

0.179148656 

03 

0.440000006 

03 

46 

0. 16549896E 

03 

0.450000006 

03 

47 

0. 14951714E 

03 

0.460000006 

03 

48 

0.13098114E 

03 

0.470000006 

03 

49 

0. 10972241E 

03 

0.48O000006 

03 

50 

0.855909586 

02 

0.490000006 

03 

51 

0.58487442E 

02 

0. 50000000E 

03 

52 

0.28366730E 

02 

0.510000006 

03 

53 

0.0 


0.518609626 

03 


WTE: Similar data are printed for each angle 

but have been omitted in this figure. 


Figure 10.- Continued. 
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12/01/69 


SULAR CELL PERFURMANCE 
RUN 2 


sample problem 2 « ARRAY CALCULATION 
flXEF) CELL characteristics 


AREA= 2.00 VREF=594.0 IKEF=60.00 ATX= 0.000350 BRS= 2.00 RS= O.ll 


ENVIRONMENTAL CELL PERFORMANCE - (COMPLETE ARRAY! 


AU = 0.5000 TC 360.0000 

NIS = ^.00 VIS = 5^,00 SCL = 0.50000000E 00 


THETAS 

THETAS 

FTHETAS 

FTHETAS 


5.00 

15.00 

55.00 

65.00 

2.67 

2.67 

2.67 

2.67 


I 

CURRENT(AMPS) 

1 

0. 14981200E 

02 

2 

0. 14981203E 

02 

3 

0.14981274E 

02 

4 

0. 14980529E 

02 

5 

0. 14980104E 

02 

6 

0. 14978915E 

02 

7 

0. 14977575E 

02 

8 

0. 149757 16E 

02 

9 

0.14972413E 

02 

10 

0. 14967646E 

02 

11 

0. 14960480E 

02 

12 

0. 14949732E 

02 

13 

0. 14933014E 

02 

14 

0. 14908407E 

02 

15 

0.14870512E 

02 

16 

0.14813605E 

02 

17 

0.14727546E 

02 

18 

0. 14597622E 

02 

19 

0.14402163E 

02 

20 

0. 14109879E 

02 

21 

0.13675360E 

02 

22 

0. 13036803E 

02 

23 

0. 12109967E 

02 

24 

0.10788912E 

02 

25 

0.89494381E 

01 

26 

0.64611082E 

01 

27 

0.34832869E 

01 

28 

0.764I9097E 

00 


25.00 35.00 45.00 

75.00 85.00 

2.67 2.67 2.67 

2.67 2.67 


VULTAGE(VOLTS) 


0.0 

O.IOUOOOOOE 01 
0.20000000E 01 
0.30000000E 01 
0.40U00000E 01 
0.50000000E 01 
0.60000000E 01 
0.70000000E 01 
0.80000000E 01 
0.90000000E 01 
O.IOOOOOOOE 02 
O.llOOOOOOE 02 
0.12000000E 02 
0.13000000E 02 
0.14UOOOOOE 02 
0.15000000E 02 
0.16000000E 02 
0.17000000E 02 
0.18000000E 02 
0.19000000E 02 
0.20000000E 02 
0.21000000E 02 
0.22000000E 02 
0.23000000E 02 
0.24U00000E 02 
0.25UO0000E 02 
0.26000000E 02 
0.27000000E 02 


KPMAX) = 0.13369101E 02 V(PMAX) = 0.20535751E 02 P(MAX) = 0.27459644E 03 


Figure 10.- Concluded 



SAMPLE PROBLEM 3 - FILTERED SOLAR CELL 


The reduction in electrical output created by a blue-red filter is sought for the cell of Sample 
Problem 1 . Table 1 lists the relative spectral response of the cell as well as the transmission of the 
coverglass filter. To illustrate the input procedure for response and transmission values that are a 
function of wavelength only, the temperature and angular effects have been neglected. 

When a value for IREF and a value of -1 for XJLTR are entered, the computer will calculate 
a value for CONST and then integrate the product of spectral response, transmission, and solar 
spectrum to obtain a true value of XJLTR for use in determining the light-generated current. To 
enable a comparison to be made with Sample Problem 1, AU = 0.50 and T = 360° K. 

The input data form, data deck listing, and computer output are shown in figures 11,12, and 13. 


TABLE 1.- TRANSMISSION AND RELATIVE RESPONSE VERSUS WAVELENGTH 


WAVE, M 

RESP 

TRANS 

WAVE, n 

RESP 

TRANS 

0.30 

0.00 

0.00 

0.85 

0.98 

0.95 

.35 

.05 

.01 

.90 

.88 

.94 

.40 

.10 

.82 

.95 

.63 

.93 

.45 

.23 

.88 

1.00 

.37 

.89 

.55 

.52 

.90 

1.05 

.14 

.90 

.65 

.80 

.87 

1.10 

.05 

.55 

.75 

.94 

.95 

1.15 

.00 

.10 

.80 

.99 

.95 

1.30 

.00 

.02 
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INPUT DATA - SOLAR CELL PERFORMANCE PROGRAM 


FORMAT INPUT DATA 
A6 RUN 5 

12A6 TITLE CARD Sample Problem 3 - Filtered Solar Cell 

6110 NAU 1 NTHETA 1 KPLOT 0 KOUT 1 INT 1 INR 1 INTD 0 INRD_0 

8F10.4 AU 0.50 

TEMP 360. 

A 1.43 

VCL 0.21 

AREA ^ TREF 305. VREF 594. IREF 60.0 XJLTR -1. ATX .00055 
BRS 2.0 RS 0.11 

THETA 0 . 

FTHETA 1 . 

NIS 0. VIS SCL 

PLOT CAR D Not lised - delete 

TABLES 

TRANSMISSION (TRANS) Table 1 


RESPONSE (RESP) Table 1 


TRANSMISSION DAMAGE (DTRANS) Not used - delete 


RESPONSE DAMAGE (DRESP) Not used - delete 


COMMENTS Filtered cell calculations illustrating single variable tables 


Figure 11.- Input data form - Sample Problem 3. 
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KUN 3 


1 

sample 

1 

PROBLEM 3 
0 

“ filtered 
1 

SOLAR CELL 
1 

1 

0 

2. 

0 

0.5 

360. 

1.43 

0.21 

2.0 

303. 

594. 

60.0 

-1. 

.00035 

.11 

0. 

1. 

0. 

16 

.30 

.35 

.40 

.45 

.55 

.65 

.75 

.80 

.85 

.90 

.95 

1.00 

1.05 

1.10 

1.15 

1.30 

.00 

.01 

.82 

.88 

.90 

.87 

.95 

.95 

.95 

.94 

.93 

.89 

.90 

.55 

.10 

.02 

16 

.30 

.35 

.40 

.45 

.55 

.65 

.75 

.80 

.85 

.90 

.95 

1.00 

1.05 

' 1. 10 

1.15 

1.30 

.00 

.05 

.10 

.23 

.52 

.80 

.94 

.99 

.98 

.88 

.63 

.37 

.14 

.05 

.00 

.00 


Figure 12.- Data deck - San 5 )le Problem 3. 


INPUT DATA RUN 3 


sample problem 3 - FILTERED SOLAR CELL 


NAU = 

1 


ntheta= 

1 


KPLOT= 

0 


KOUT = 

1 


AREA = 

0.20000000E 

01 

TREF = 

0.30300000E 

03 

VREF = 

0.59400000E 

03 

IREFs 

0.60000000E 

02 

XJLTRs-O. lOOOOOOOE 

01 

ATX = 

0.34999987E-03 

BRS* 

0.200000006 

01 

RS« 

0.10999995E 

00 


ABSOLUTE RESPONSE/RELATIVE RESPONSE = 0.53334 


Figure 13.- Output 


Sample Problem 3. 



12/01/69 


SOLAR CELL performance 
RUN 3 


SAMPLE PROBLEM 3 - FILTERED SOLAR CELL 
FIXED CELL characteristics 


AREA= 2.00 VR£F=S9A.O 1REF=60.00 ATX= 0.000350 BRS= 2.00 RS= O.U 


ENVIRONMENTAL CELL PERFORMANCE - (SINGLE CELL) 


AU= 0.50 THETA = 0.0 TC =360.0 JLTR = 0.19393A80E 00 VCL =0.21 

IL = 0.2215A022E 03 10 = 0. 2083 1837E-02 VUC = 0.5139233<»E 03 A = 1.43 


I CURRENT(MA) VOLTAGE(MV) 


1 

0.22154022E 

03 

0.0 


2 

0.22154022E 

03 

0.1 0000000 E 

02 

3 

0.22154022E 

03 

0.20000000E 

02 

4 

0.22154022E 

03 

0.30000000E 

02 

5 

0.22154022E 

03 

0.400000006 

02 

6 

0.22154022E 

03 

0.50000000E 

02 

7 

0.22152837E 

03 

0.600000006 

02 

8 

0.22152837E 

03 

0.700000006 

02 

9 

0.22152837E 

03 

0.800000006 

02 

10 

0.22151492E 

03 

0.900000006 

02 

11 

0.22150801E 

03 

0. lOOOOOOOE 

03 

12 

0.22150705E 

03 

O.llOOOOOOE 

03 

13 

0.22148849E 

03 

0.12000000E 

03 

14 

0.22147491E 

03 

0.13000000E 

03 

15 

0.22146600E 

03 

0. 140000006 

03 

16 

0.22143657E 

03 

0.150000006 

03 

17 

0.22140990E 

03 

0.16000000E 

03 

18 

0.22138596E 

03 

0.170000006 

03 

19 

0.22133458E 

03 

0. 18000000E 

03 

20 

0.22128215E 

03 

0.19000000E 

03 

21 

0.22121649E 

03 

0.20000000E 

03 

22 

0.22113751E 

03 

0.2 10000006 

03 

23 

0.22103134E 

03 

0.220000006 

03 

24 

0.22089728E 

03 

0.23000000E 

03 

25 

0.22073517E 

03 

0.24000000E 

03 

26 

0.22054456E 

03 

0.250000006 

03 

27 

0.22028679E 

03 

0.260000006 

03 

28 

0.21996173E 

03 

0.270000006 

03 

29 

0.21956877E 

03 

0.280000006 

03 

30 

0.21908218E 

03 

0.290000006 

03 

31 

0.21846545E 

03 

0.300000006 

03 

32 

0.21769563E 

03 

0.310000006 

03 

33 

0.21673543E 

03 

0. 320000006 

03 


rt ->i litirtocc 

03 


03 

35 

0.21405034E 

03 

0.34000000E 

03 

36 

0.21220084E 

03 

0.350000006 

03 

37 

0.209907536 

03 

0.360000006 

03 

38 

0.20707069E 

03 

0.370000006 

03 

39 

0.203571826 

03 

0.380000006 

03 

40 

0. 199271156 

03 

0.390000006 

03 

41 

0.194007146 

03 

0.400000006 

03 

42 

0.187595726 

03 

0.410000006 

03 

43 

0.179831286 

03 

0.420000006 

03 

44 

0.170490976 

03 

0.43000000E 

03 

45 

0.159338936 

03 

0.440000006 

03 

46 

0.14613478E 

03 

0.450000006 

03 

47 

0.130644336 

03 

0.460000006 

03 

48 

0.112649126 

03 

0.470000006 

03 

49 

0.919502416 

02 

0.480000006 

03 

50 

0.68411240E 

02 

0.490000006 

03 

51 

0.419166116 

02 

0.500000006 

03 

52 

0.124122406 

02 

0.510000006 

03 

53 

0.0 


0.513923346 

03 


I(PMAX) = 0.19800998E 03 V(PMAX) = 0.39260571E 03 P(MAX) = 0.77741623E 02 

Figure 13.- Concluded. 
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SAMPLE PROBLEM 4 - SIMULATED SOLAR ENERGY 


Laboratory tests are scheduled for the cell described in Sample Problem 3. The spectral energy 
output of the solar simulator is tabulated in table 2. What differences in cell performance can be 
expected under simulated conditions for 0.50 AU? 

Information on the spectral energy distribution of the source is placed in the first three data 
statements of subroutine RESPE. These contain a list of wavelengths, the ratios of energy in the 
bandwidths centered at these wavelengths to the total energy output, and the total number of 
wavelengths listed. Figure 14 illustrates the placement of the statements. 

The input data form, data deck listing, and computer output are shown in figures 15,16, and 1 7. 


TABLE 2.- SOLAR SIMULATOR SPECTRUM 


WAVE^ 

LSOLAR^ 

WAVE® 

lsolar’’ 

0.31 

0.0174 

0.83 

0.0181 

.33 

.0056 

.85 

.0283 

.35 

.0072 

.87 

.0066 

31 

.0183 

.89 

.0084 

.39 

.0243 

.91 

.0139 

.41 

.0286 

.93 

.0113 

.43 

.0298 

.95 

.0080 

.45 

.0320 

.97 

.0090 

.47 

.0339 

.99 

.0086 

.49 

.0358 

1.01 

.0093 

.51 

.0342 

1.03 

.0071 

.53 

.0311 

1.05 

.0044 

.55 

.0306 

1.07 

.0061 

.57 

.0312 

1.09 

.0061 

.59 

.0296 

1.11 

.0089 

.61 

.0284 

1.13 

.0121 

.63 

.0278 

1.15 

.0068 

.65 

.0264 

1.17 

.0065 

.67 

.0251 

1.19 

.0139 

.69 

.0249 

1.21 

.0091 

.71 

.0252 

1.23 

.0052 

.73 

.0218 

1.25 

.0046 

.75 

.0213 

1.27 

.0046 

.77 

.0213 

1.29 

.0075 

.79 

.81 

.0218 

.0177 

1.31 

.0039 


^Center wavelength of 0.02 n band. 

^Energy in 0.02 ju band/total energy output. 
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C SUBROUTINE RESPE - INTERPOLATION 

SUBROUTINE RESPE I THE T t TC f AUC f FJL TR r K T ) 
REAL IREF,LSnLAR,LREF 


DIHENS ION XYZ (3) ,XYZ1(3) , X YZ 2 ( 3 ) , X Y Z3 ( 3 ) ,NXYZ (3 ) ,NXYZ1 (3) , 

1 NXYZ2(3),NXYZ3|3),TINV(150),T0INV(150),RINV(150), 

2 RDINV(150),TARGa260l,T0ARG(lZ60l,RARG(1260I tR0ARG(126OI, 

3 PAROERO ) ,HAVE (90) tLSOLAR (90 ) , RWAVE ( 90 ) t LRE F ( 90 ) t 

A TRANS(90),DTRANS(90),RESP(90) »DRESP(90) ,FACT0R(90) 

COMMON /COMl/AREA, IREF,TREF, INT,INRtINTO,lNRO 
C SUBSTITUTION OF XENON - KRYPTON SPECTRUM FOR SULAR SPECTRUM 
C A CENTER WAVELENGTH OF 0.02 MICRON BANDWIDTH UF ENERGY SOURCE 

DATA WA VE/0,3 1,0. 33 1 0.35,0. 37,0. 39, O.A I t0.^t3T0.^5,O.A7,O.A9,0. 51 1 

1 0.53,0.55,0.57,0.59,0.61,0.63,0.65,0.67,0.69,0.71,0.73, 

2 0.75,0.77,0.79,0.81,0.83,0.85,0.87,0.89,0.91,0.93,0.95, 

3 0.97,0.99,1.01,1.03,1.05,1.07,1.09,1.11,1.13,1.15,1.17, 

A 1.19,1.21,1.23,1.25,1.27,1.29,1.30/ 

C FRACTION OF TOTAL ENERGY IN BANDWIDTH OF ENERGY SOURCE 
C (RELATIVE DISTRIBUTION) 

DATA LSOLAR/. 017^, .0056, .007 2,. 0183, .02A3, .0286, .0298, .0320, .0339, 

1 .0350, . 0342, .0311, .0306, .0312, .0296, .0284 ,. 0278 ,. 0264 , 

2 .0251,.0249,.0252,.0218,.0213,.0213,.0218,.0177,.0iai, 

3 .0283,.0066,.0084,.0139,.0113,.0080,.0090,.0086,.0093, 

4 .0071 , .0044, .0061, .0061 , . 0089 ,. 0 12 1 ,. 0068 , .0065 , .0 1 39 , 

5 .0091 , .0052, .0046, .0046,.0075, .0039/ 

C NW IS THE NUMBER OF VDLUES IN THE WAVE AND LSOLAR ARRAYS AND MUST 
C BE redefined IF THE ARRAY SIZE IS CHANGED 

DATA NW/51/ 

C center wavelength of 0.02 MICRON BANDWIDTH UF ENERGY SUURCE USED TO 
C OBTAIN IREF 

DATA RWAVE /O. 31, 0.33, 0.35,0. 37, 0.39,0. 4 1,0. 43, 0.45,(3.47,0. 49,0.51, 


C 

C 


C 

C 


1 0.53,0.55,0.57,0.59,0.61,0.63,0.65,0.67,0.69,0.71,0.73, 

2 0.75,0.77,0.79,0.81,0.83,0.85,0.87,0.89,0.91,0.93,0.95, 

3 0.97,0.99,1.02,1.06,1.10,1.14,1.18,1.22,1.26,1.30/ 
fraction of total ENERGY IN BANDWIDTH OF ENERGY SOURCE 

(RELATIVE DISTRIBUTION) 


DATA 

LREF /.0103 

,.0154 

, .0167 

,.0182 

.0174 

.0267 

.0270, .0310 

.0310 

1 

.0290, 

0280, 

0270, 

.0280, 

.0270, 

02 70, 

.0250, .0240, 

0240, 

2 

.0220, .0210, 

.0200, 

.0190, 

.0180, 

0180, 

.0160, .0160, 

0160, 

3 

.0140, 

0140, 

.0140, , 

.0130, 

,0120,, 

.0110, , 

.0110, .0110,, 

.0210, 

4 

.0180, 

0170, 

.0160, , 

.0150, 

.0140,, 

.0130, 

.0120/ 



NP IS THE NUMBER OF VALUES IN THE RWAVE AND LREF ARRAYS AND MUST 
BE redefined if THE ARRAY SIZE IS CHANGED 


data NP/43/ 

t TIRR is total IRRADIANCE IN MV/SO. CM. OF SOURCE, LREF, USED TO 

C OBTAIN IREF 


Figure 14.- 


Energy spectrum substitution in subroutine RESPE. 
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INPUT DATA - SOLAR CELL PERFORMANCE PROGRAM 


FORMAT INPUT DATA 
A6 RUN 4 

12A6 TITLE CARD Sample Problem 4 - Simulated Solar Energy 

6110 NAU I NTHETA I KPLOT 0 KOUT I INT 1 INR I INTD 0 INRD_0 
8F10.4 AU 0.50 

temp 360. 

A 1.45 

VCL 0.21 . 

AREA 2 _^ TREF 505. VREF 594. I-REF 60.0 XJLTR -I. ATX .00035 
BRS 2.0 RS O.ll 

theta 0.0 


FTHETA 1.0 

NIS 0 . VIS SCL 

PLOT CAR D Not used - delete 

TABLES 

TRANSMISSION (TRANS) Table 1 
RESPONSE (RESP) Table 1 


TRANSMISSION DAMAGE (DTRANS) Not lised - delete 


RESPONSE DAMAGE (DRESP) Not used - delete 


COMMENTS Illustrates input procedure for energy sources other than the sun 


Figure 15.- Input data form - Sample Problem 4. 
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RUN 4 



SAMPLE 

PRObLEM ^ - 

SIMULATED 

1 

0.5 
360. 
1 .^3 
0.21 

1 

0 

1 

2.0 

0* 

1. 

0. 

lb 

303. 

59<^. 

60.0 

.30 

.35 

.40 

.45 

.85 

.90 

.95 

1.00 

.00 

.01 

.82 

.86 

,95 

16 

.9^ 

.93 

.89 

.30 

.35 

.40- 

.45 

.85 

.90 

.95 

1.00 

.00 

.05 

.10 

.23 

.98 

.88 

.63 

.37 


SOLAR tNERGY 

I I 


-1. 

.00035 

2. 

.11 

.55 

.65 

.75 

.80 

I.U5 

I.IO 

1.15 

1.30 

.90 

.87 

.95 

.95 

.90 

.55 

.10 

.02 

.55 

.65 

.75 

.80 

1.05 

1.10 

1.15 

1.30 

.52 

.80 

.94 

.99 

.14 

.05 

.00 

.00 


Figure 16.- Data deck - Sample Problem 4. 


INPUT DATA RUN A SAMPLE PRUBLtM A - SIMULATED SULAR ENERGY 


NAO = 

1 


MTHETA= 

1 


KPLOT= 

0 


KUUT = 

1 


AREA = 

0.200000UOE 

01 

TKEF = 

0.3O3OUUO0E 

03 

VKEF = 

0.5940UOUOE 

03 

IREF = 

0.6000UO00E 

02 

XJLTR = ‘ 

-0. lOOOOOOOE 

01 

ATX = 

0.34999987E-03 

bKS = 

0.2000UOO0E 

01 

RS = 

0. 10999995E 

00 


absolute response/relative response = 0.5333A 


Figure 17.- Sample Problem 4. 
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l^f/03/69 


SULAR CeiL PEKFURMANCt 
RON 


SAMPLE PRUhLEM ^ ' SIMULATED SULAR ENERGY 
FIXED CELL characteristics 


AREA= 2,00 VREF=59^.0 IKEF=60.00 ATXs 0.0003SO BRS= 2.00 RS* 0.11 


ENYIRONMeNTAL CELL PERFORMANCE - (SINGLE CELL! 


AU= 0.50 THETa = 0.0 TC =360.0 JLTR « 0.2165I097E 00 VCL =0.21 

IL = 0.2A732941E 03 10 = 0, 20586273E-02 VUC = 0.5193391 IE 03 A = 1.A3 


CURRENT (MAI 


VULTAGE(MV1 


1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

lA 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2 ^ 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3A 

35 

36 

37 

38 

39 
AO 
A1 

42 

43 

44 

45 

46 

47 

48 

49 

50 

52 

53 


0.24732941E 03 
0.24732941E 03 
0.247329416 03 
0.24732941E 03 
0.247329416 03 
0.24732941E 03 
0.24731679E 03 
0.24731679E 03 
0.24731679E 03 
0.2473U263E 03 
0.24729533E 03 
0.247294166 03 
0.24727478E 03 
0.24726048E 03 
0.24725079E 03 
0.247220U9E 03 
0.24719199E 03 
0.24716603E 03 
0.24711266E 03 
0.24705743E 03 
0.246908276 03 
0.246904306 03 
0.24679326E 03 
0.246654316 03 
0.24648752E 03 
0.246275076 03 
0.24601642E 03 
0.245676456 03 
0.245254706 03 
0.244740706 03 
0.244091636 03 
0.243281376 03 
0.242271036 03 
0.241012506 03 
0.239447026 03 
0.237502936 03 
0.235093586 03 
0.232115076 03 
0.228444096 03 
0.22393616E 03 
0.218424126 03 
0.211719016 03 
0.20361093E 03 
0.19387337E 03 
0.182268336 03 
0.168556756 03 
0.15250624E- 03 
0.133903926 03 
0.11255859E 03 
0.S8345016E 02 
0.61157120E 02 
U.30951050E 02 
0.0 


0.0 

0.100000006 02 
0.200000006 02 
0.300000006 02 
0.400000006 02 
0.500000006 02 
0.600000006 02 
0.700000006 02 
0. 800000006 02 


.900000006 02 
0.100000006 03 
0. 110000006 03 
0. 120000006 03 
0.130000006 03 
0.140000006 03 
0.150000006 03 
0.160000006 03 
0.170000006 03 
0. 180000006 03 
0.190000006 03 
0.2U0000006 03 
0. 210000006 03 
O.22OOO0006 03 
0.230000006 03 
0.240000006 03 
0.25000000E 03 
0.260000006 03 
0.270000006 03 
0.28000000E 03 
0.290000006 03 
0.300000006 03 
0.310000006 03 
0.32U000006 03 
0. 330000006 03 
0.340000006 03 
0. 350000006 03 
0.36000000E 03 
0. 370000006 03 
0.380000006 03 
0. 390000006 03 
0. 400000006 03 
0.410000006 03 
0.42O0UO00E 03 
0.430000006 03 
0. 440000006 03 
0.450U00006 03 
O.46OOO00OE 03 
0.470000006 03 
0.480000006 03 
0. 490000006 03 
0. 500000006 03 
0.510000006 03 
0.519339116 03 


I(PMAX) = 0.220997B5E 03 


V(PMAX) = 0.395597416 03 


P(MAX) = 0.874282386 02 


Figure 17.- 


Concluded 


SAMPLE PROBLEM 5 - COMPLETE MISSION ANALYSIS 


The output characteristics of the solar array described in Sample Problem 2 are to be 
determined for a typical heliocentric orbit with perihelion at 0.4 AU. The cells are similar to those 
of Sample Problem 3 and have transmission and response values as shown in tables 3 and 4. To 
illustrate the bivariant use of the tables, the relative response is shown to be a variable of wavelength 
and incident angle. Transmission is a function of wavelength, incidence angle, and temperature to 
demonstrate the trivariant capability of the program. 

The reduced response of the cells at high angles of incidence is due to an experimentally 
determined decrease in output due to factors other than the cosine reduction in incident energy and 
the spectral transmission shift of the filters. 

Data on transmission damage to coverglasses and response damage to the solar cells have been 
compiled into tables 5 and 6 for a typical heliocentric orbit. Note that it is the difference between 
the degraded and undegraded values that is used in the tables. 

The value of CONST has previously been determined to be 0.533 and may be entered 
directly by setting IREF equal to -1 and entering 0.533 for XJLTR. 

Calculations will be made at 1.0, 0.8, 0.6, and 0.4 AU. The corresponding panel temperatures 
are expected to be 285°, 320°, 367°, and 450° K, respectively. Values for A and VCL are obtained 
from figures 1 and 1 8, respectively. Neither plotting nor detailed output is desired. 

The input data form, data deck listing, and computer output are shown in figures 1 9, 20, and 21 . 
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TABLE 3 .- FILTER TRANSMISSION ( TRANS ) VERSUS WAVELENGTH ( KWAVE ), 
INCIDENCE ANGLE ( KTHET ), AND TEMPERATURE ( KTEMP ) 




200 ° - 340 ° 

K 





500 ° 

K 



X.M 

0 ° 

15 ° 

30 ° 

45 ° 

60 ° 

90 ° 

0 ° 

15 ° 

30 ° 

45 ° 

60 ° 

m 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


.01 

.01 

.01 

.01 

.01 

.10 

.01 

.01 

.01 

.10 

.01 

.15 


.02 

.02 

.02 

.09 

.20 

.40 

.05 

.06 

.12 

.40 

.60 


380 

.06 

.06 

.12 

.40 

.60 

.85 

.40 

.42 

.57 

.75 

.80 

.85 

390 

.42 

.42 

.57 

.75 

.85 

.80 

.83 

.80 

.82 

.84 

.88 

.90 

400 

.82 

.80 

.82 

.84 

.82 

.75 

.82 

.80 

.82 

.84 

.82 

.75 

410 

.85 

.85 

.85 

.84 

.80 

.73 

.83 

.83 

.83 

.81 

.78 

.71 

460 

.88 

.87 

.86 

.84 

.80 

.75 

.86 

.85 

.84 

.82 

.78 

.73 

550 

.90 

.89 

.88 

.86 

.82 

.77 

.88 

.87 

.86 

.84 

.80 

.75 

620 

.91 

.90 

.88 

.80 

.76 

.70 

.89 

.88 

.85 

.78 

.74 

.68 

640 

.90 

.89 

.83 

.86 

.80 

.76 

.88 

.87 

.81 

.84 

.78 

.74 

650 

.87 

.86 

.84 

.87 

.82 

.78 

.85 

.84 

.82 

.85 

.80 

.76 

670 

.92 

.91 

.91 

.90 

.85 

.80 

.89 

.87 

.87 

.86 

.82 

.78 

700 

.92 

.91 

.93 

.92 

.87 

.85 

.90 

.89 

.91 

.90 

.85 

.83 

750 

.95 

.94 

.92 

.92 

.90 

.87 

.92 

.91 

.89 

.89 

.87 

.84 

850 

.95 

.94 

.95 

.90 

.85 

.77 

.92 

.91 

.92 

.87 

.82 

.74 

950 

.93 

.93 

.88 

.85 . 

.70 

.50 

.90 

.90 

.85 

.82 

.67 

.48 

1000 

.89 

.88 

.90 

.56 

.25 

.05 

.86 

.85 

.77 

.53 

.10 

.02 

1050 

.90 

.91 

.60 

.15 

.02 

.02 

.40 

.39 

.37 

.23 

.02 

.02 

1100 

.55 

.10 

.10 

.04 

.03 

.04 

.02 

.02 

.02 

.02 

.02 

.02 

1150 

.10 

.04 

.03 

.04 

.08 

.10 

.02 

.04 

.06 

.08 

.10 

.12 

1200 

.02 

.02 

.02 

.05 

.10 

.20 

.02 

.05 

.10 

.15 

.25 

.30 

1250 

.02 

.02 

.03 

.10 

.30 

.50 

.02 

.10 

.20 

.40 

.65 

.70 

1300 

.02 

.03 

.05 

.25 

.65 

.70 

.05 

.30 

.40 

.65 

.70 

.75 
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TABLE 4.- RELATIVE SPECTRAL RESPONSE (RESP) VERSUS WAVELENGTH (MWAVE) 

AND INCIDENCE ANGLE (MTHETA) 


X.M 

0° 

45° 

60° 

75° 

85° 

90° 

0.30 

0 

0 

0 

0 

0 

0 

.35 

.05 

.049 

.048 

.044 

.038 

.033 

.40 

.100 

.099 

.096 

.088 

.075 

.065 

.45 

.230 

.228 

.221 

.202 

.173 

.150 

.55 

.520 

.514 

.499 

.458 

.390 

.338 

.65 

.80 

.792 

.768 

.704 

.600 

.520 

.75 

.940 

.931 

.902 

.827 

.705 

.611 

.80 

.990 

.980 

.950 

.871 

.743 

.644 

.85 

.980 

.970 

.941 

.862 

.735 

.637 

.90 

.880 

.871 

.845 

.774 

.660 

.572 

.95 

.630 

.624 

.605 

.554 

.473 

.410 

1.00 

.370 

.366 

.355 

.326 

.278 

.241 

1.05 

.140 

.139 

.134 

.123 

.105 

.091 

1.10 

.050 

.049 

.048 

.044 

.038 

.033 

1.15 

.000 

.000 

.000 

.000 

.000 

.000 

1.30 

.000 

.000 

.000 

.000 

.000 

.000 
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TABLE 6.- RESPONSE DAMAGE VERSUS WAVELENGTH (JWAVE) AND AU (JAU) 

[0 = 0° to 90° ] 


AU 


Damage 

1.0 

0.3 

0.000 


.8 

.000 


1.1 

.000 


1.3 

.000 

0.6 

.3 

.002 


.8 

.01 


1.1 

.002 


1.3 

.002 

0.4 

.3 

.002 


.8 

.020 


1.1 

.006 


1.3 

.004 

0.2 

.3 

.006 


.8 

.030 


1.1 

.010 


1.3 

.007 




VCL, mV/mA 


.40 



Ql I I I I I I 1 1 

■ 0 4 8 12 16 20 24 28 

Solar constants 


Figure 18,- Solar constants. 


35 



INPUT DATA - SOLAR CELL PERFORMANCE PROGRAM 


FORMAT INPUT DATA 
A6 RUN 5 


12A6 TITLE 

CARD 

Sample Problem 5 

- Complete Mission Analysis 


6110 NAU_4 

NTHETA 9 KPLOT 0 

KOUT 0 

INT 5 

INR 2 INTO 2 

INRD 2 

8F10.4 AU 

1.0 

0.8 0.6 

0.4 




TEMP 

285. 

320. 367. 

450. 




A 

2.0 

1.90 1.62 

1.22 




VCL 

0.30 

0.28 0.26 

0.24 




AREA 

2.0 TREF 303. VREF 594 

I REF 

-1. 

XJLTR 0.533 ATX 

.00035 


BRS 2.0 RS 0.11 





THETA 

5. 

15. 25. 

35. 

45. 

55. 65. 

75. 


85. 






FTHETA 

2.667 

2.667 

2.667 2.667 

2.667 

2.667 

2.667 2.667 

2.667 

NIS 

4. 

VIS 54. SCL 

0.5 




PLOT CARD_ 

Not used - delete 






TABLES 


TRANSMISSION (TRANS) Table 3 


RESPONSE (RESP) Table 4 


TRANSMISSION DAMAGE (DTRANS) Table 5 


RESPONSE DAMAGE (DRESP) Table 6 


COMMENTS Array calculation including radiation damage and direct input 
of CONST 


Figure 19.- Input Data Form - Sample Problem 5. 
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RUN 5 


sample problem 5 “ complete mission analysis 


4 

9 

0 

0 

3 

2 

2 

2 

1.0 

0.8 

0.6 

0.4 





285 . 

320 . 

367 . 

450 . 





2.0 

1 .9 

1.62 

1.22 





0.30 

0.28 

0.26 

0.24 





2.0 

303 . 

594 . 

- 1 . 

0.533 

0,00035 

2.0 

0.11 

5 . 

15 . 

25 . 

35 . 

45 . 

55 . 

65 . 

75 . 

85 . 








2.667 

2.667 

2.667 

2.667 

2.667 

2.667 

2.667 

2.667 

2.667 









54 . 

0.5 






24 

6 

3 






.300 

,350 

.370 

.380 

.390 

.400 

.410 

.460 

.550 

.620 

.640 

.650 

.670 

.700 

.750 

.850 

.950 

l.OOO 

1.050 

1.100 

1.150 

1.200 

1.250 

1.300 

0 . 

15 . . 

30 . 

45 . 

60 . 

90 . 

200 . 

340 , 

500 . 








.0 

.01 

.02 

.06 

.42 

.82 

.85 

.88 

.90 

.91 

.90 

.87 

.92 

.92 

.95 

.95 

.93 

.89 

.90 

.55 

.10 

.02 

.02 

.02 

.0 

.01 

.02 

.06 

.42 

.80 

.85 

.87 

.89 

.90 

.89 

.86 

.91 

.91 

.94 

.94 

.93 . 

.88 

.91 

.10 

.04 

.02 

.02 

.03 

.0 

.01 

.02 

.12 

.57 

.82 

.85 

.86 

.88 

.86 

.83 

.84 

.91 

.93 

.92 

.95 

.86 

.90 

.60 

.10 

.03 

.02 

.03 

.05 

.0 

.01 

.09 

.40 

.75 

.84 

.84 

.84 

.66 

• 80 

.86 

.87 

.90 

.92 

.92 

.90 

.85 

.56 

.15 

.04 

.04 

.05 

.10 

.25 

.0 

.01 

.20 

.60 

.85 

.82 

.80 

.80 

.82 

.76 

.80 

.62 

.85 

.87 

.90 

.85 

.70 

.25 

.02 

.03 

.08 

.10 

.30 

.65 

.0 

.10 

.40 

.85 

.80 

.75 

.73 

.75 

.77 

.70 

.76 

.78 

.80 

.85 

.87 

.77 

.50 

.05 

.02 

.04 

.10 

.20 

.50 

.70 

.0 

.01 

.02 

.06 

.42 

.82 

.85 

.88 

.90 

.91 

.90 

.87 

.92 

.92 

.95 

.95 

.93 

.89 

.90 

.55 

.10 

.02 

.02 

.02 

.0 

.01 

.02 

.06 

.42 

.80 

.85 

.87 

.89 

.90 

.89 

.86 

.91 

.91 

.94 

.94 

.93 

.68 

.91 

.10 

.04 

.02 

.02 

.03 

.0 

.01 

.02 

.12 

.57 

.82 

.85 

.86 

.88 

.88 

.63 

.84 

.91 

,93 

.92 

.95 

.88 

.90 

.60 

.10 

.03 

.02 

.03 

.05 

.0 

• 01 

.09 

.40 

.75 

.84 

.64 

.84 

.86 

.60 

.86 

.87 

.90 

.92 

.92 

.90 

.85 

.56 

.15 

.04 

.04 

.05 

.10 

.25 

.0 

.01 

.20 

• 60 

.85 

.82 

• 60 

.80 

.82 

.76 

.60 

.82 

.85 

.87 

.90 

.85 

.70 

.25 

.02 

.03 

• 08 

;io 

• 30 

.65 

.0 

.10 

.40 

.85 

.80 

.75 

.73 

.75 

.77 

.70 

.76 

.78 

.80 

,85 

.87 

.77 

.50 

.05 

.02 

.04 

.10 

.20 

.50 

.70 

.0 

.01 

.05 

.40 

.83 

.82 

.83 

.86 

.86 

.89 

.88 

.85 

.89 

.90 

.92 

.92 

.90 

.86 

.40 

.02 

.02 

• 02 

.02 

.05 

.0 

.01 

.06 

.42 

.80 

.80 

.83 

,85 

.87 

.88 

.87 

.84 

.87 

.89 

.91 

.91 

.90 

.85 

.39 

.02 

.04 

.05 

• 10 

.20 

.0 

.01 

.12 

.57 

.82 

.82 

.83 

.64 

.86 

.85 

• 61 

.82 

.87 

.91 

.89 

.92 

.85 

,77 

.37 

.02 

.06 

. 10 

.20 

.40 

.0 

.01 

.40 

.75 

.84 

.84 

.61 

.82 

.84 

.78 

.84 

.85 

.66 

,90 

.89 

.87 

.82 

.53 

.23 

.02 

.08 

.15 

.40 

.65 

.0 

.01 

• 60 

.80 

.88 

.82 

.78 

.78 

.80 

.74 

.78 

.80 

.82 

.85 

.87 

.82 

.67 

.10 

.02 

.02 

.10 

.25 

.65 

.70 

.0 

.15 

• 80 

.65 

.90 

.75 

.71 

.73 

.75 

.68 

.74 

.76 

.78 

.83 

.84 

.74 

.48 

.02 

.02 

.02 

.12 

.30 

.70 

.75 

16 

6 

2 






.30 

.35 

• 40 

.45 

.55 

.65 

.75 

.80 

.85 

.90 

.95 

l.OO 

1.05 

1.10 

1.15 

1.30 

0 . 

45 . 

60 . 

75 . 

65 . 

90 . 

100 , 

500 . 

.000 

.050 

.100 

.230 

.520 

.800 

.940 

.990 

.980 

.880 

.630 

.370 

.140 

.050 

.000 

.000 

.000 

.049 

.099 

.228 

.514 

.792 

.931 

.980 

.970 

.871 

.624 

.366 

.139 

.049 

.000 

.000 

.000 

.04 8 

.096 

.221 

.499 

• 768 

.902 

.950 

.941 

.845 

.605 

.355 

.134 

.048 

.000 

.000 

.000 

.044 

.088 

.20 2 

.458 

.704 

.827 

.871 

.862 

.774 

.554 

.326 

.123 

.044 

.000 

.000 

.000 

.038 

.075 

.173 

.390 

.600 

.705 

.743 

.735 

.660 

.473 

.278 

.105 

.038 

.000 

.000 

.000 

.033 

.065 

.150 

.338 

.520 

.611 

.644 

.63 7 

.572 

.410 

.241 

.091 

.033 

.000 

.000 

.000 

.050 

.100 

.230 

.520 

.800 

.940 

.990 

.980 

.880 

.630 

.370 

.140 

.050 

.000 

.000 

.000 

.049 

.099 

.228 

.514 

.792 

.931 

.980 


Figure 

20.- 

Data deck - 

San^le 

Problem 

5. 
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► 970 

.871 

.624 

.366 

.139 

.049 

.000 

.000 

► 000 

.048 

.096 

.221 

.499 

.768 

.902 

.950 
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,1)00 
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.000 

.000 

► 000 

.038 

.075 

.173 

.390 

.600 

.705 

.743 

► 735 

.660 

.473 

.278 

.105 

.038 

.000 

.000 

,000 

.033 

.065 

.150 

.338 

.520 

.611 

.644 

► 637 

.572 

.410 

.241 

.091 

.033 

.000 

.000 

8 

2 
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1. 
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,0 

90. 

.2 

.4 
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.004 
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.0 

,005 

.005 

.016 

.009 

.014 

.004 

.0 

.0 

,008 

.008 

.024 

.015 

.021 

.005 

.0 

.0 

,008 

.008 

.024 

.015 

.021 

.005 

.0 

.0 

,02 

.018 

.048 

.029 

,040 

.012 

.001 

.001 

,02 

.018 

.048 

.029 

.040 

.012 

.001 

.001 

,12 

.11 

.39 

.23 

.345 

.085 

.015 

,01 

,12 

.11 

.39 

.23 

.345 

.085 

.015 

.01 
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1.3 

.0 
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.4 

,6 

1. 







,0 

.0 

.0 

.0 
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.0 

.002 

.01 
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.01 
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.02 

.006 

.004 

.004 
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.004 
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.007 

.006 
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Figure 20.- Concluded. 


input UATA run 5 


SAMPLE PRUbLtM 5 - COMPLETE MISSION ANALYSIS 


NAU= ^ 

NTHETA= 9 
KPLOT= 0 
KUUT= 0 

AREA= 0.20000000E Ol 
TKEF= 0.30300UUOE 03 
VKEF= O.59A0U0OOE 03 
lkEF= 0.0 

XJLTR= 0.53299999E 00 
aTX= 0.3^999987E-03 
bRS= 0.2000UOUOE Ol 
RS= 0.10999995E 00 


AbSOLUTE RESPOnSE/RELATIVE KESPUNSE = 0.53300 


Figure 21.- Output - Sample Problem 5, 



10/21/70 


SOLAR CPLL PERFORMANCE 
RUN 5 

SAMPLE PROBLEM 5 - COMPLETE MISSION ANALYSTS 

fixed cell characteristics 

APFA=» 2.00 VRFF=5O^,0 TREF=-1 ,00 ATX= 0,000350 RR S= 2,00 RS= 0,11 
ENVIRONMENTAL CELL PERFORMANCE - (COMPLFTF ARRAVf 


AU = 1,0000 

TC = 

285,0000 


NTS = ^,00 

VIS 

= SA.OO 

SCL = O.SOOOOnOOE 


THETAS » 

5.00 

1 5,00 

25, 00 

35,00 

45,00 

THETAS = 

55 ,00 

6F ,00 

75,00 

85.00 


FTHETAS 

?,F7 

2,57 

7.57 

7.57 

2.67 

FTHFTAS - 

7.F7 

2.57 

2.57 

2.57 



1 

CUPRENTIAMPSI 

VOLTAOEI VOLTS! 

1 

0, 25333708E 

01 

0.0 


2 

0.26333213E 

01 

0, lOOOOOOOE 

01 

3 

0.2533321 8E 

01 

n.2oooooooF 

01 

4 

0.26333227F 

01 

0. 30000000E 

01 

5 

0.253T3277F 

01 

0.40000000F 

01 

6 

o. 25333208E 

01 

o.soonnoooE 

01 

7 

0.253331 89E 

01 

o.6oonooooF 

01 

8 

0.25333094E 

01 

0.70000O00F 

01 

9 

0. 26332941 E 

01 

0. 80000000E 

01 

10 

0.25332584E 

01 

0.90000000F 

01 

1 1 

0.26332216E 

01 

O.lOOOOOnOF 

02 

17 

0.26331405F 

01 

0. iinoooooF 

02 

13 

0-76325895E 

01 

0. 12000O00F 

02 

14 

O. 26325750E 

01 

0, 13O0OO00E 

02 

1 5 

0.2632O333F 

01 

0, 14000000F 

02 

16 

0.26313410E 

01 

0. 15On0000F 

02 

17 

0.25303644E 

01 

0. 150000O0F 

02 

18 

0.26289206E 

01 

0.17000000E 

02 

19 

0, 26258940E 

01 

0. I8000000E 

02 

20 

0.26238489F 

01 

0. 1900nr>00F 

0? 

21 

0.2619387RE 

01 

0.20000000E 

0? 

22 

0.26129723E 

01 

0. 21 OOOOOOF 

02 

23 

P.25035835E 

01 

0.22000000E 

02 

74 

0, 25999705F 

01 

0.73000000E 

02 

25 

0.25702677E 

01 

0. 24000000F 

02 

26 

0,2541 7750E 

01 

O.25000O00E 

07 

77 

0. 24996090E 

01 

0. 75000000E 

0? 

28 

0.24388313E 

01 

n,27000000F 

0? 

29 

0.23506231E 

01 

O.28000000E 

02 

30 

0.22227736F 

01 

0. 29000000F 

02 

31 

0.20380P17F 

01 

0.300000 OOF 

02 

32 

0. 17725830F 

01 

0,31000000E 

02 

33 

0.13924150E 

01 

0. 32000O00F 

02 

34 

0.85285771 E 

00 

0,33000n00E 

07 

35 

0.13250685E 

00 

0.34000000F 

02 

0. 740274S7E 01 

VTPMAXI = 0.27464981E 

02 PIMAXl 

= 0.65998932E 0? 


Figure 21.- 


Continued 



10 / 21 /70 


SrHAR CELL PEPPHRMANCE 
RUN •> 


SAWPLF PROPLFM 5 - COH»LFTF MISSION ANALYSIS 
FIYEO CELL CHARACTFRISTTCS 


APEA= 2,00 VRFF=594.n rffCF=-I,00 ATX= ^^,0(^0?50 0RS= 2,00 RS= 0,11 


FNVTRONMFNTAL cell performance - (COMPLETF ARRAY) 


All = O.ROnO TC = ■5?0.0000 


NIS = 4.00 VIS = *>4.00 Sr.L = O.SOOOOOOOE 00 


thetas = 

5.00 

15.00 

25.00 

35.00 

45.00 

THETAS = 

55,00 

65.00 

75.00 

85.00 


FTHETAS = 

2.67 

2.67 

2.67 

2.67 

2.67 

FTHETAS = 

2.67 

2.67 

7.67 

2.67 



! 

CURRENTl AMPS > 

VPLTAGFl VOLTS) 

1 

0.46R466OSE 

01 

0.0 


'7 

0.468466R5F 

01 

O.IOOOOOOOF 

01 

3 

0.46R46704E 

01 

0. 200000006 

01 

4 

0.46346695F 

01 

0,300000006 

01 

5 

0.46846695E 

01 

0,400000096 

01 

6 

0.46846342E 

01 

0.500000006 

01 

7 

0.45841526E 

01 

0.60000000F 

01 

8 

0.468384176 

01 

0.700000006 

01 

0 

0.468331056 

01 

0.80000000F 

01 

I 0 

0.46824973F 

01 

0.90000000F 

01 

11 

0.468147326 

01 

O.IOOOOOOOF 

02 

12 

0.467995936 

01 

O.llOOOOOOF 

0? 

13 

0.467778596 

01 

0.17000000E 

02 

14 

0.46746655F 

01 

0. 13000000F 

02 

15 

0. 46704655 F 

01 

0.14000000F 

02 

16 

0.46640615E 

01 

0. 15000000F 

0? 

17 

0.465532406 

01 

0.160000006 

0? 

1 8 

0.46425228F 

01 

0,170000006 

02 

19 

0.46248856F 

01 

0, 18000000E 

02 

70 

0.4590 0738F 

01 

0.190000006 

02 

21 

0.45631227F 

01 

0.2OOOO0006 

02 

72 

0.45116863F 

01 

0.210000006 

02 

73 

0.443902976 

01 

0.22000000F 

02 

74 

0.43T61O69F 

01 

0.230000006 

02 

75 

0.419057226 

01 

0,24000000F 

02 

76 

0.39S51522E 

OJ 

0,250000006 

02 

77 

0.369649226 

01 

0.26000000F 

02 

?fl 

0.329294116 

01 

0,270000006 

02 

?o 

0. 27335272F 

01 

0. 280000006 

02 

30 

0.196420576 

01 

0. 29000000F 

02 

31 

0.921952136 

00 

0,300000006 

02 


I(pMAX) = n.419‘?70R?E 01 VIPMAX) = 0.23970291F 0? P(MAX) = 0.l0n57326E 0? 

NOTE: Similar data for 0.6 and 0.4 AU are printed but 

have been eliminated in this figure. 


Figure 21.- Concluded. 
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APPENDIX B 


DESCRIPTION OF DATA CARD INPUT 



Control Cards for 360 j 61 Version 
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Control Cards for 7094 Version 
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MAIN PROGRAM 


m 

Variable 

name 

Columns 

Format 

Description of variables 

1 

R(l) 

1-6 

A 6 

Run number 


R(l), R(2) 

1-8 

2A4^ 


2 

TITLE 

1-72 

12A6(18A4)^ 

Run identification 

3 

NAU 

1-10 

8110 

Number of variables in AU, T, A, 





and VCL arrays 


NTHETA 

11-20 


Number of variables in THETA and 





FTHETA arrays 


KPLOT 

21-30 


1 = plot (S-C 4020) 
0 = omit plot 


KOUT 

31-40 


1 = detailed single-cell output 
0 = no single-cell output 


INT 

41-50 


Input code for transmission table 





= 0 Omit 

= 1 Singularly variant (use TAINT) 
= 2 Bivariant (use TVIN) 

= 3 Trivariant (use TVIN) 


INR 

51-60 


Input code for response table 





(same as INT) 


INTD 

61-70 


Input code for transmission 





damage table (same as INT) 


INRD 

71-80 


Input code for response damage 





table (same as INT) 

4 b 

AU 

1 - 10 , etc. 
1 - 10 , etc. 

8F10.4 

Solar distance, AU 

5b 

T 

8F10.4 

Temperature, °K 

6 b 

A 

l- 10 ,etc. 

8F10.4 

Curve factor used in diode equation 





(eq. ( 6 )) 

7 b 

VCL 

l- 10 ,etc. 

8F10.4 

Current coefficient used in open 





circuit voltage equation (eq. 5)), 
mV/mA 

8 

AREA 

1-10 

8F10.4 

Single-cell surface area, cm^ 


TREE 

11-20 


Reference temperature, °K 


VREF 

21-30 


Reference voltage, mV 


IREF 

3M0 


Reference current (AMZ), mA 


^Format for IBM 360/67. 

^Eight values per card; number of cards depends on NAU. 
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Card 

Variable 

Columns 

Format 

Description of variable 

number 

name 



8 

XJLTR 

41-50 


If (INT + INR) = 0 





0 Conversion constant used in 
equation (1) in A/mW 
= 0 Compute XJLTR from 





equation (2) 





If (INT + INR) f 0 





f -1 Set CONST = XJLTR and solve 
equation (3) 

= -1 Calculate CONST from 





equation (4) and then solve 
equation (3) for XJLTR 


ATX 

51-60 


Temperature coefficient in current 





density equation (eq. 2)), 
Al/I per °K 


BRS 

61-70 


Temperature coefficient in open 





circuit voltage equation (eq. (5)), 
mV/°K 


RS 

71-80 


Cell series resistance, fZ 

9C 

THETA 

1-10, etc. 

8F10.4 

Solar incidence angles, deg 


FTHETA 

1-10, etc. 

8F10.4 

Number of strings at angular 





orientation THETA 

11 

NIS 

1-10 

8F10.4 

Number of cells in a parallel group 


VIS 

11-20 


Number of parallel groups wired in 





series to form a string 


SCL 

21-30 


String circuit loss, V 

12^ 

PLOT 

14-20 


Date 


INFORMATION 

25-32 


Name 



34-37 


Installation 



40-42 


Building number 

SUB 

ROUTINE RESPE 

Transmission ta 

ble (required only if INT f 0) 

13 

KWAVE 

1-10 

3110 

Number of wavelength variables 


KTHET 

1 1-20 


Number of incidence angle variables 


KTEMP 

21-30 


Number of temperature variables 

14® 

XTAB 

1-10 

8F10.4 

Array of wavelengths, angles, and 





temperatures. Array size is 
KWAVE + KTHET + KTEMP 

15® 

WTAB 

1-10 

8F10.4 

Transmission values corresponding 



1 l-20,etc. 


to XTAB. Array size is 
KWAVE*KTHET*KTEMP 


^Eight values per card; number of cards depends on NTHETA. 

^Only if KPLOT = 1 . 

^Number of cards depends on array size. See description of Subroutine AL TVIN in appendix E. 


45 








Card 

Variable 

Columns 

Format 

Description of variable 

number 

name 



Response tables (required only if INR ^ 0) 

16 

MWAVE 

1-10 

3110 

Number of wavelength variables 


MTHET 

11-20 


Number of incidence angle variables 


MTEMP 

21-30 


Number of temperature variables 

17® 

XTAB2 

1-10, etc. 

8F10.4 

Array of wavelengths, angles, and 





temperatures. Array size is 
MWAVE + MTHET + MTEMP 

18® 

WTAB2 

1-10, etc. 

8F10.4 

Response values corresponding to 





XTAB2. Array size is 
MWAVE*MTHET*MTEMP 

Transmission damage tables (required only if INTD 4 0) 

19 

LWAVE 

1-10 

3110 

Number of wavelength variables 


LTHET 

11-20 


Number of incidence angle variables 


LAU 

21-30 


Number of AU variables 

20® 

XTABl 

l-10,etc. 

8F10.4 

Array of wavelengths, angles, and 





AU’s. Array size is 
LWAVE + LTHET + LAU 

21® 

WTABl 

l-10,etc. 

8F10.4 

Transmission damage values 






corresponding to XTAB 1 . Array 
size is LWAVE* LTHET* LAU 

Response damage tables (required only if INRD 4 0) 

22 

JWAVE 

1-10 


Number of wavelength variables 


JTHET 

11-20 


Number of incidence angle variables 


JAU 

21-30 


Number of AU variables 

23® 

XTAB3 

l-10,etc. 


Array of wavelengths, angles, and 





AU’s. Array size is 
JWAVE + JTHET + JAU 

24® 

WTAB3 

l-10,etc. 


Response damage values 





corresponding to XTAB3. Array 
size is JWAVE*JTHET*JAU 


^Number of cards depends on array size. See description of Subroutine AL TVIN in appendix E. 
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APPENDIX C 


Variable 

A 

AA 

ABC ( 
ABSC J 

AREA 

ATX 

AU 

AUC 

BL 

BU 

BRS 

DENOM 

DMON 1 \ 
DMON 2 1 
DMON 3 I 
DMON 4 I 
DMON 5 I 
DMON 6 / 



EO 

ET 

FMAX 


PROGRAM SYMBOL LIST 

Main Program 

Description 

curve factor used in equations (6) and (7) 
single value in A array 

plotting titles for the X-axis 

single solar cell surface area 

temperature coefficient in equation (1) 

solar distance in astronomical units 

single value in the AU array 

lower bound on X in subroutine AL ROOT 

upper bound on X in subroutine AL ROOT 

i 

temperature coefficient in equation (5) 
denominator of equation (6) 

variables used by AL TAINT which must be set to zero initially and whenever 
the argument array is altered 

numbers that specify the degree of precision expected of the root 
determination (AL ROOT) 

equals El 

equals E2 

interpolation variable used in subroutine PMAX 


47 



Variable 


Description 


FTHETA 

FUNC 

GS 

II 

I 

ILOAD 

IL 

INR 

INRD 

INT 

INTD 

10 

IREF 

JCDREF 

JCD 

KINT 

KOUT 

KPLOT 

KSV 

KT 

LAST 

NAU 

NER 


number of strings at each orientation (eq. (10)) 
product of current and voltage in subroutine PMAX 
initial guess (AL ROOT) 

number of voltage and current pairs used to describe I-V curve 
see ILOAD, SCURR 

current at a given voltage (equivalent to I in equations (6) and (8)) 

light generated current 

input code for response tables 

input code for response damage tables 

input code for transmission tables 

input code for transmission damage tables 

reverse saturation current (eq. (6)) 

reference short circuit current (eq. (2)) 

factor used in computing XJLTR 

short circuit current density (eq. (1)) 

sum of table input codes used to regulate use of subroutine RESPE 
input code regulating detailed output 
input code regulating plotting 
counter used in array calculations 

constant used by subroutine RESPE to determine if tables should be read 
variable to test end of data 
size of AU array 
error code (AL TAINT) 
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Variable 

NIS 

NRR 

NSV 

NTHETA 

ORD 1 ' 
ORD 2 
ORD 3 ‘ 
ORD 4, 

PILOAD 

PI 

PSP 

PVOLTM 

PWRMAX 

R 

RS 

SCL 

SCURR 

SFUNC 

SI 

SMAX 

SVOLT 

T 

TC 

THETA 

TITLE 


Description f 

number of cells connected in parallel (eq. (8)) 
error code (AL TAINT, AL ROOT) 
counter used in array calculations 
size of NTHETA and FTHETA arrays 

plotting titles for y-axis 

ILOAD*(-1.0) 

7T, 3.1415927 

specific power (power/area) 
not used 
maximum power 
array containing title for each run 
series resistance (eq. (6)) 

'sfnhg'cifcuiribss(eqr(9)) — 

string current variables (equivalent to I in eq. (10)) 

string current x voltage 

variable used in string current calculations 

current at maximum power 

string voltages (equivalent to V in eq. (11 )) 

solar cell operating temperature 

single variable in T array 

solar incidence angles in degrees 

title for printer output 



Variable 


Description 


TREF 

V 

VCL 

VIS 

VOC 

VOLTM 

VREF 

XJLTR 

XL 

XMAX 

XU 

YL 

YU 


reference temperature (eq. ( 1 )) 

see VOLTM, SVOLT 

current coefficient (eq. (5)) 

number of cells connected in series (eq. (9)) 

open circuit voltage (eq. (5)) 

voltage variable used in obtaining I-V curve from equation (6) (equivalent to 
V in eqs. (6) and (9)) 

reference voltage (eq. (5)) 

conversion constant (eq. (1)) 

lower limit of x-axis for plotting 

voltage at maximum power 

upper limit of x-axis for plotting 

lower limit of y-axis for plotting 

upper limit of y-axis for plotting 


CONST 

DMON 
DMON 1 
DMON 2 
DMON 3 ) 

DRESP 

DTRANS 

FACTOR 

lERR 

JAU 


Subroutine RESPE 
proportionality factor in equation (3) 

variable used by AL TAINT which must be set to zero initially and whenever 
the argument array is altered 

response damage corresponding to a specific wavelength (eq. (3)) 

transmission damage corresponding to a specific wavelength (eq. (3)) 

variable used in calculation of XJLTR 

error code from (AL TVIN) 

number of solar variables in response damage tables 
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Variable 


Description 


JTHET 

JWAVE 

K1 

K2 

K3 

K4 

KTEMP 
KTHET 
KWAVE 
LAU 
LREF 
LSOLAR 
LTHET 
-L.WA_YE_ 
MTEMP 
MTHET 
MWAVE 
NER 
NP 
NW 
NW 1 
NW 2 
NW 3 


number of incidence angle variables in response damage tables 
number of wavelength variables in response damage tables 
number of variables in transmission array (TINV) 
number of variables in transmission damage array (TDINV) 
number of variables in response array (RINV) 
number of variables in response damage array (RDINV) 
number of temperature variables in transmission tables 
number of incidence angle variables in transmission tables 
number of wavelength variables in transmission tables 
number of solar distance variables in transmission damage tables 
source spectral energy distribution used in equation (4) 
solar spectral energy distribution (L\) used in equation (3) 
number of incidence angle variables in transmission damage tables 
number of w avelen gth variables in transmission damage tables 
number of temperature variables in response tables 
number of incidence angle variables in response tables 
number of wavelength variables in response tables 
error code (AL TAINT) 

number of variables in the RWAVE and LREF data statements 
number of variables in the WAVE and LSOLAR data statements 
number of variables in the transmission array (TARG) 
number of variables in the transmission damage array (TDARG) 
number of variables in the response array (RARG) 



Variable 
NW 4 
NXYZ 

NXYZ 1 

NXYZ 2 

NXYZ 3 

RARG 

RDARG 

RDINV 

RESP 

RINV 

RWAVE 

SUM 

TARG 

TDARG 

TDINV 

TINV 

TIRR 

TRANS 

WAVE 


Description 

number of variables in the response damage array (RDARG) 

array containing the number of wavelengths, angles, and temperatures 

in TINV table 

array containing the number of wavelengths, angles, and temperatures 

in TDINV table 

array containing the number of wavelengths, angles, and temperatures 

in RINV table 

array containing the number of wavelengths, angles, and temperatures 

in RDINV table 

array of response values corresponding to RINV table 

array of response damage values corresponding to RDINV table 

array of wavelengths or wavelengths, angles, and temperatures for response 
damage 

response variable corresponding to a wavelength as a result of table look-up 

array of wavelengths or wavelengths, angles, and temperatures for response 

array of wavelengths corresponding to LREF array (eq. (4)) 

variable used in computing XJLTR 

array of transmission values corresponding to TINV table 

array of transmission damage values corresponding to TDINV table for 
transmission 

array of wavelengths or wavelengths, angles, and temperatures for transmission 
damage 

array of wavelengths or wavelengths, angles, and temperatures for transmission 
constant used to compute CONST 
transmission variable 

array of wavelengths corresponding to LSOLAR (eq. (3)) 
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Variable 


Description 


XYZ array of wavelength, angle, and temperature arguments for table look-up for 

transmission 

XYZ 1 array of wavelength, angle, and temperature arguments for table look-up for 

transmission damage 

XYZ 2 array of wavelength, angle, and temperature arguments for table look-up for 

response 

XYZ 3 array of wavelength, angle, and temperature arguments for 'table look-up for 

response damage 
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2 


0=>.N THETA' 
k (No. ongles) . 
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SUBROUTINE PMAX 





APPENDIX E 


PROGRAMMER’S INFORMATION 


The program is coded in FORTRAN IV in two versions. One version is designed to run on an 
IBM 7094 computer with optional graphical output on an SC-4020 plotter. The decks were 
compiled and the program was validated on the IBM 7094 computer under the IBJOB processor 
using the IBSYS operating system, version 13. The second version is designed to run on an 
IBM 360/67 computer, with optional graphical output on an SC-4020 plotter. These decks were 
compiled on an IBM 360/67, using Release 17, FORTRAN H with OPT = 2 and the program was 
validated on an IBM 360/67 computer using OS, Release 17. 

The program should operate with only minor changes in control cards in either compiler. The 
amount of core available on the IBM 7094 computer was 32,768 words. On the IBM 360/67 
computer 1 58,000 bytes of core was available. 

If a SC-4020 plotter is not available, the IBM 7094 version of the program may be altered as 
follows: 

1 . Omit these cards 
EXTERNAL SYSACC 
IF(KPLOT. EQ. 0)GO TO 2 
CALL READIN (SYSACC) 

SO^TALLEOFTV " “ 

2. Change 
CALL EXIT 
to 

50 CALL EXIT 

3. Omit data card supplied for CALL READIN 

4. Substitute a dummy subroutine for AL AICR 

The IBM 360/67 version of the program should be altered in the following manner: 



1 . Omit these cards 
IF(KPLOT. EQ. O)G0 TO 2 
CALL READIN (SYSACC) 

2. Change 

5000 CALL EQ)FTV 
to 

5000 STOP 

3. Omit data card supplied for 
CALL READIN 

4. Substitute a dummy subroutine for AL AICR 
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KPLOT = 0 SKIP PLOTTING 
KPLOT = 1 PLOT J 

KOUT = 0 SKIP DETAILED OUTPUT FOR CURRENT VS VOLTAGE (THETA) 
KOUT = 1 PRINT DETAILED OUTPUT FOR CURRENT VS VOLTAGE (THETA) 
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SUBROUTINE PM AX ( PWRMAX tXMAX ,FMAX ) 
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200 FORMAT ( 1H0,12HTAIMT ERR0R=,I2) 
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SOLAR CELL PERFORMANCE PROGRAM - KIRKPATRICK 
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213 FORMAT (IH ,25X, IHI ,9X , 13HCURRENT ( AMPS ) ,9X, 14HV0LTAGE (VOLTS) // ) 
DO 21 1=1, NSV 



21 CONTINUE 
36 CONTINUE 

CALL TAINT ( S VOLT , SCURR , XMAX , SMAX ,NSV , 2 t NER , DMON 1 ) 

WRITE (6,205) SMAX , XMAX , PWRMAX 

205 FORMATdHO/ ,1H0, 5X,9HI(PMAX) = , El 5 . 8 , 5X , 9HV ( PMAX ) = , E 15 . 8 , 5X , 8HP I 
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APPENDIX G 


DESCRIPTION OF LIBRARY SUBROUTINES 
IDENTIFICATION 


AL TVIN, trivariate table look-up 
FORTRAN IV, MAP-coded 

Ames modification of SHARE library routine SC TRIV 


PURPOSE 


This subroutine is used to evaluate the function W = f(X,Y,Z), and the partial 
derivatives 9W/9X, 9W/9Y, 9W/9Z, by linear interpolation, where W is given as a tabulated 
function of X, Y, and Z. 


USAGE 


This subroutine is called by use of the statement 
CALL TVIN(W, XYZ, XTAB, WTAB, lERR, NXYZ, PARDER) 
where 

W contains the result of the table look-up 

XYZ is an array that contains the arguments for table look-up 

XYZ(l) contains X 
XYZ(2) contains Y 
XYZ(3) contains Z 

XTAB is an array that contains the tables of X, Y, and Z values upon which the search is to 
be made; 

XTAB(l) - XTAB(NX) contains the values for Xj , Xj, . . ., Xj^^ 

XTAB(NX + 1 ) - XTAB(NX + NY) contains the values for Yj , Yj, . . ., Yjsjy 
XTAB(NX -I- NY -I- 1) - XTAB(NX + NY + NZ) contains the values 
for Zi, Z 2 , . . ., Zn 2 
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WTAB 


is an array that contains functional values: 


WTAB( 1 ) - WTAB(NX) contains W(X, , Yj ,Zj ) - WCXj^x- Yi ,Zj ) 

WTAB(NX + 1) -WTAB(2*NX) contains W(Xi ,Yj,Z, )- W(X nx,Y 2 ,Z, ) and so 
forth, with Z varying last, until 
WTAB(NX*NY*NZ) contains W(Xnx.YnY>^NZ) 

lERR is an error return code such that 

lERR = 1 , normal return 
= 2, X is too small 
= 3, X is too large 
= 4, Y is too small 
= 5, Y is too large 
= 6, Z is too small 
= 7, Z is too large 

NXYZ is an array that contains the number of variables in XTAB: 

NXYZ(l) contains NX 
NXYZ(2) contains NY 
NXYZ(3) contains NZ 

PARDER is an array that contains the partial derivatives: 

PARDER(l) contains aW/3X 
PARDER(2) contains 9W/9Y 
PARDER(3) contains 9W/9Z 


RESTRICTION 


The values entered in the three tables in XTAB must increase monotonically. 


COMMENTS 


This subroutine performs the same function as the subroutine AL TRIV, which is also a 
modification of SC TRIV. The present version is more convenient to use, however, as tables are 
loaded in forward order rather than reverse order. 


METHOD 


The following interpolation formulas are used: 
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The quantity X^. is Xj or Xj+j, whichever is closer to X. If X is halfway between Xj and Xj+, , 
then Xj+i is used as X^.. If X = XTAB(NX), for example, the last entry in the table, 
then Xj+i is XTAB(NX) and Xj is XTAB(NX-l). The above discussion applies for Y^. and 
Zq also. 


IDENTIFICATION 


AL TAIN(TAINT), table look-up and interpolation 
FORTRAN IV, MAP-CODED 
Written by V. L. Sorensen 


PURPOSE 


This subroutine is used to evaluate one or more functions of an argument X, when the 
functions and the argument are given in tabular form. 

USAGE 

The CALL statement for the IBM 7040/7094 DCS for this subroutine is 
CALL TAINT(XTAB, FTAB, X, FX, N, K, NER, DMON, GTAB, GX, HTAB, HX, . . .) 
and the CALL statement for the IBM 360/67 is 
CALL TAINT(XTAB, FTAB, X, FX, N, K, NER, DM0N) 
where 
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XTAB 


is the name of the array that contains the values of the argument of the functions, 
which must increase or decrease monotonically 


are the names of arrays that contain the values of the functions of the argument 
in XTAB 


is the value of the argument for which function values are sought 


are the values of the functions as obtained by the subroutine 


is the number of tabulated values, and so the arrays XTAB, FTAB, . . . , are of 
dimensions at least N;N>K 

K is the order of interpolation; K < 9 

NER is an error code: 

NER = 1, normal return 
NER = 2, K>9orN<K 

NER = 3, adjacent values within the array XTAB are equal 

DMON is a variable that must be set to zero initially, and whenever XTAB is altered. 

Each XTAB within a given program should be assigned a unique DMON symbol. 
The contents of DMON are altered by the subroutine. 


RESTRICTIONS 


The table XTAB must be monotonic, and adjacent values may not be equal. The 
order K must be smaller than or equal to 9, and N must be greater than K. 


METHOD 


A binary search is used to find the region of the table where the desired result is to be found. 
Aitken’s method of interpolation is then used to determine the precise functional values 
corresponding to the given value of the independent variable. Aitken’s method is a computational 
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scheme for evaluating the Lagrangian interpolation polynomial without having to compute 
polynomial coefficients. The development of the Lagrangian interpolation method will be omitted 
here, as details may be found in many standard texts (i.e., see ref. 1). The Aitken method (see 
ref. 2) consists in the repetitive use of the formula 


X) - - X) 

m+n^ n-i 


The superscripts and subscripts are indices, and never exponents, in this formula. For a given value 
of the order K, the above formula is applied for n in the range 1, 2, . . ., K, and, for each n, the 
values for m are taken in the order 0, 1, . . ., K-1. The restriction n + m>K applies. Thus, for 
linear interpolation (K = 1 ), this process is carried out with a single application of the above general 
formula, for example, 

yJ(X^ - X) - yJcXq - X) 

" xT^x^ 

For the case K = 3, the process is repeated six times. The six iterations give 


1 Ygcx, - X) - y;cXo - X) 
^ ■ Xi - Xo 

1 yJCX2 - X) - Y°(Xq - X) 

'^2 ' *2 - Xo 

„1 xjcx, - X) - y"cX(, - X) 

'^3 xT^ov 

„2 y|(X 2 - X) - Y^(Xi - X) 

■ x5“n<i 

2 yJcXs - x) - y‘cx, - X) 


„3 X^CXj - X) - Y^CXj - X) 

X3 ~ Y ., _ 


n = 1, m = 0 


n = 1 , m = 1 


n = 1, m = 2 


n = 2, m = 0 


n = 2, m = 1 


n = 3, m = 0 


Note that the fourth and fifth iterations depend on the first, second, and third, and that the sixth 
and final one depends on the fourth and fifth. It is evident that the Aitken method may be 
characterized as the repetitive application of a linear interpolation process. 
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IDENTIFICATION 


AL SC42, S-C 4020 microfilm recorder plotting package 
FORTRAN IV 360/50 

Written by North American Aviation, obtained from UAIDE 


PURPOSE 


AL SC42 is a package of FORTRAN IV subroutines used in preparing input tapes for 
the S-C 4020 Microfilm Recorder (Plotter) operated by Lockheed Missiles and Space Company, 
Sunnyvale, California. The S-C 4020 produces plotted and printed output on vellum and/or 
microfilm. The use of this machine offers great advantages for jobs having large volumes of plotted 
output. 


USAGE 


An excellent manual for users of the S-C 4020 is supplied by Lockheed (ref. 1), and copies 
are available from the Computation Division. Usage described in this manual is for an IBM 7094 
computer; however, only a few differences exist for usage with an IBM 360/50. These differences 
are listed at the end of this write-up. 

The Ames user must include the following in his program: 

1. CALL READIN * 

This statement must appear prior to calling any of the ALSC42 subprograms. He must supply a 
data card to be read by READIN, containing the following information: 


Columns 

Item 


1-68 

Name of user, mail stop, and any 



other alphanumeric information 

* 

69-72 

Job ID 

♦ 

73-80 

Date 

* 


The information contained on this card is used by ALSC42 for the leading and trailing 
identification frames. The information in columns 69 through 80 is also used for 
the ID information in the upper right-hand corner of each frame. The CALL READIN statement 
also generates the leading identification frame. 
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2. CALL E0FTV 


This statement is used to generate the trailing ID frame on the plot tape followed by an end of file, 
and also writes on tape 6 the number of frames plotted. 

3. Necessary control cards 

INCLUDE DECKS(ALSC42) * 

//G0.FTO1OO1 DD UNIT=TAPE, V0LUME=SER=ASSIGN, X * 

// DCB=(DEN=2, TRTCH=C, RECFM=U, BLKSIZE=400), X * 

// DUSP=(NEW,KEEP), LABEL=(1,BLP) X * 

The storage required by ALSC42 is approximately 45,000 bytes. * 

The *’s that appear at the right of the page indicate differences between this write-up and the 7094 
write-up of ALSC42. 


CURRENT MACHINE TIME RATES 


As of this date, machine usage charges for the S-C 4020 Microfilm Plotter are as follows: 

(1 ) $90.00 per elapsed machine hour time plus 

(2) (a) 5 cents per frame microfilm (1 6 mm) 

(b) 25 cents per frame hard copy (vellum) 

(c) 27 cents per frame for combination of hard copy and microfilm 

(3) $25.00 minimum charge per job submission shall apply to combined cost elements of (1) 
and (2) above. Generally, a job submission will include plot tapes from several users. 


REFERENCES 


North American Aviation, Inc., S-C 4020 Manual. 

Lockheed Missiles and Space Co., Sunnyvale, California, December 1965. 


ALSC42 PLOTTING PACKAGE FOR 360/50 


General Description 

The S-C 4020 manual that is distributed describes subroutine usage for the IBM 7094 version 
of ALSC42. The IBM 360/50 version of ALSC42 provides the user with essentially all the 
capability described in this manual. The differences between the two versions are listed below. 
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1 . The following routines are no longer available: 

CHAINV 

CLEANV 

DSHLNV, see item 21 

ERRLNV 

GRACS 

IDFRM, IDFRM., IDFRMS 

SCCTAB 

SCP00L 

SETFRV 

VLAGX 

VLAGY 

2. No routine may have an optional calling sequence of different length than the normal 
calling sequence. When this restriction appeared to decrease capability, it was circumvented by 
adding a new routine using the optional calling sequence (see PRINTV and TYPEV). 

3. No routine may pass locations of variables instead of values (see SCERRV, SERSAV, 
SERREV). No capability is apparently lost by this restriction, although the usage of the three 
routines mentioned is changed somewhat. 

4. APLOTV - The characters to be plotted (stored in the array MARKET) are assumed to be 
represented by integers (that are right adjusted) if NC is positive, and by characters (that are left 
adjusted) if NC is negative. Thus, if APLOTV is called with a Hollerith argument for MARKET, 
or if MARKET consists of characters read in with an A1 field, NC must be negative. 

5. BNBCDV - The results returned by BNBCDV occupy two words, thus BCDWD must be 
an array dimensioned at least two. The six characters are right adjusted. 

6. E0FTV - No longer rewinds and unloads the plot tape. 

7. FRAMEV - Must always be called with an argument. 

8. NOFRV - May be called with one argument only 

CALL NOFRV(NOFRM). 

Returns the current frame count number. 

9. POINTV - May be called with three arguments only 

CALL POINTV(X,Y, ±NS). 

10. PRINTV - May be called with four arguments only 

I CALL PRINTV(N,BCDTXT,IX,IY). 
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The “type current point” feature is utilized by a new routine TYPEV that is the same 
as PRINTV with two arguments. 

1 1 . RESETV - May not have an argument. Use FRAMEV if comer marks and frame counting 
are desired. 

1 2. SCERRV - The arguments represent values of the scaling error indicators instead of 
locations. 


CALL SCERRV(KX,KY). 

Stores the current values of these indicators in KX and KY. Thus it is necessary to 
call SCERRV after using NXV and/or NYV instead of before as described in the manual. See 
also SERSAV and SERREV below. 

13. SERSAV and SERREV - These retrieve and restore the values of the scaling error 
indicators, respectively. Thus, SCERRV and SERSAV both perform the same function 
while SERREV allows the programmer to set the indicators to whatever value he desires. 

14. TABLrV-Only three tables of vector characters are implemented at the present 
time-TABLlV, TABL2V, and TABL3V. The structure of vector character tables is somewhat 
different than described in the manual in that each 12-bit pattern occupies one-half word (16 bits) 
of the core. 

1 5. XAXISV and YAXISV - May be called with three arguments only 

CALL XAXISV(IX,IY,NSTPT). 

16. TYPEV - New routine to utilize “type current point” feature of the S-C4020. Usage is 
identical to the description in the manual of PRINTV with two arguments 

CALL TYPEV(N,BCDTXT). 

1 7. SHFTIV - New routine to perform logical one register, off-the-end shift 

CALL SHFTIV(AIN,AOUT,NS). 


where 

AIN input 
AOUT result 

NS number of bit positions to shift, left if plus, right if minus 
Vacated bit positions are filled with zeros. 
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1 8. SHFT2V - New routine to perform logical two register, off-the-end shift 

CALL SHFT2V(AIN,BIN,AOUT,BOUT,NS) 

Arguments similar to SHFTIV. On a left shift (NS+), bits are shifted through bit position 1 
of BIN into bit position 32 of AIN and vice versa for right shift. Vacated positions are filled with 
zeros. 

1 9. ANDV - New routine to perform logical product of two quantities. 

CALL ANDV(A,B,Q 

performs C = A * B. 

20. ORAV - New routine to perform logical turn of two quantities. 

CALL ORAV(A,B,C) 

performs C = A + B. 

21. DOTLNV, INCRV - New routine for plotting a dashed line. 

CALL D0TLNV(IX1,1Y1,IX2,IY2) 


where 

IXl.IYl start coordinates (integer values) 

1X2, 1 Y2 stop coordinates (integer values) 

The routine INCRV is used to supply the line and space size. This is set normally at 8 for the line 
and 4 for the space. These values may be altered, however, by calling INCRV prior to DOTLNV. 

CALL INCRV(IL,IS) 


where 

IL desired length of line in raster units (integer value) 

IS desired size of space in raster units (integer value) 

If IL and IS are such that the stop coordinates (IX2,IY2) terminate in the space portion of a dash 
segment, the lost line portion extends to the stop coordinates. Illegal values of IL or IS are ignored 
by INCRV. DOTLNV will then use the last used values of IL and IS. 
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IDENTIFICATION 


AL AICR, S-C 4020 Microfilm Recorder Plotting Program 
FORTRAN IV 360/50 

Written by North American Aviation, obtained from Lockheed Missiles and Space Company 


PURPOSE 


AL AICR is a general purpose subroutine for writing a magnetic tape for displaying graphical 
output on the S-C 4020. Hard copy and/or microfilm can be obtained. The subroutine will 
compute scale factors, generate the printing of the scales, generate grid lines and axes, and generate 
titles and labels from information supplied by the user. It uses the S-C 4020 plotting package 
(ALSC42) for preparing the magnetic tape that is used as input by the S-C 4020 plotter operated 
by Lockheed. This plotter offers great advantages for jobs having large volumes of output. 


USAGE 


This subroutine is entered by the use of the statement: 

CALL AICRT3(KX,KY,X,Y,NP,ND.NV,NS,NC,T,A,0,NF,NG,DCX.DCY,NX0,XL,XU,NY(i),YL,YU) 
where 

KX 0 for linear display in X, 1 for log display in X 

KY 0 for linear display in Y, 1 for log display in Y 

X name of the array containing the values of X 

Y name of the array containing the values of Y 

NP number of points in X and Y array to be plotted 

ND interval at which points will be displayed; for example, every point =1, every tenth 
point = 10 

NV 2 if one wishes to join points plotted with a straight line, 1 for points only 

NS 2 if the data are to be sorted as a function of X before plotting (to enable the connecting 

of points by vectors), 1 if data are not to be sorted 

NC decimal equivalent of the character of the point to be plotted (see figs. 1-10 attached), for 
example, NC = 38 for a circle 
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T 1 2-word alphanumeric title to be displayed at the top of the graph (see Special Features) * 

A 13-word abscissa title (see Special Features) * 

(p 1 3-word ordinate title (see Special Features) * 

NF 1 if a new frame is required, 2 if overlay is required 

NG 1 scale labeling for non-log plots is determined by the magnitude of the numbers. If 

any X meets the following condition: |X|>10® or if all X are in the range 
-10® <X< 10®, the labeling would automatically be in the E-label notation. 2 non-log 
labeling is always in the E-label notation. 


DCY* 1 coarseness of the X and Y meshes (see Special Features) 

NX0 2 override search for X-upper and X-lower, 1 let routine determine X-upper and X-lower 
(see Special Features) 

XL X-lower, if needed (NX0 = 2) 

XU X-upper, if needed (NX0 = 2) 

NY0 2 override search for Y-upper and Y-lower, 1 let routine determine Y-upper and Y-lower 
YL Y-lower, if needed (NY0 = 2) 

YU Y-upper, if needed (NY0 = 2) 

The Ames user must provide the following: 

1. CALL READIN * 


This statement must appear prior to calling the AL AICR subroutine. He must supply a data card, 
to be read by READIN, containing the following identification information: * 

Columns Item 


1-68 

Name of user, mail stop, and any 



other alphanumeric information 

* 

69-72 

Job ID 

* 

73-80 

Date 

* 


The information supplied on the card is used by the S-C 4020 for the leading and trailing 
identification frames. The information supplied in 69-80 will also appear in the upper right-hand 
comer of each frame. The CALL READIN statement generates the leading identification frame. 
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2. CALL E0FTV 


This statement is used to generate the trailing identification frame on the plot tape followed by an 
end of file and also writes on tape 6 the number of frames plotted. 

3. Necessary control cards: 

INCLUDE DECKS(ALSC42,ALAICR) * 

//G0.FTO1FOO1 DD UNIT=T APE, V<Z>LUME=SER= ASSIGN, X * 

// DCB=(DEN=2,TRTCH=C,RECFM=U,BLKSIZE=400), X * 

// DISP=(NEW,KEEP),LABEL=(1,NL) * 

The storage required by ALSC42 and AL AIRCR is approximately 5 1 ,000 bytes. * 

The *’s that appear at the right of the page indicate differences between this write-up and the 7094 
write-up of ALAICR. 


SPECIAL FEATURES 


The recommended value for .DCX and DCY is 16.0. If more grid lines are necessary, decrease 
the values of DCX and DCY. For fewer lines, increase the values of DCX and DCY. The entire 
array of X and Y being plotted must be in the core when this routine is executed. 

The feature of sorting the functions before plotting has one drawback: after the functions 
have been sorted, the original order is totally destroyed. This feature should not be used unless the 
plotting of these data is the last function of the code or unless the order of the data is not 
important. 

If (XU/XL - 1 ) <.0001 , or if XU = XL = 0, no graph will be drawn as the function is constant 
or almost constant (similarly for YU and YL). 

The titles used should be centered in the 48-character array to insure centering of the titles on 
the graph. The abscissa and ordinate title can be as large as 52 characters. 

If one uses the upper and lower limit search provided by AICRT3, care should be exercised in 
defining the variable or constant used for XL, XU, YL, or YU as it will be changed by the 
subroutine. 

NG is used to permit a floating point notation labeling of linear plots. This is provided 
because the maximum fixed point number that can be used is 999999, and the smallest greater than 
zero is .000001 , similarly to negative values. To provide meaningful labeling of data outside this 
range, the floating point notation has been made available. 

The limit of the number of cycles in log plotting is 10. Do not try to plot more than a 
lOX 10 log X log plot. 
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RESTRICTION 


This routine should not be used to plot a point at a time. 


CURRENT MACHINE TIME RATES 


As of this date, machine usage charges for the S-C 4020 Microfilm Plotter are as follows: 

(1 ) $90.00 per elapsed machine hour time plus 

(2) (a) 5 cents per frame microfilm (16 mm) 

(b) 25 cents per frame hard copy (vellum) 

(c) 27 cents per frame for combination of hard copy and microfilm 

(3) $25.00 minimum charge per job submission, shall apply to combined costs elements of (1) 
and (2) above. 


REFERENCES 


North American Aviation, Inc., S-C 4020 Manual. 

Lockheed Missiles and Space Co., Sunnyvale, California, December 1 965. 
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